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ABSTRACT 

The interaction of electromagnetic waves of wavelength X with 
periodic structures of spatial period A are studied. The emphasis of 
the work is on Bragg interactions where X 2A/N and the Bragg order 
N takes on the values 1,2,... . An extended coupled waves (ECW) 
theory is developed for the case N ^ 2 and the results of the theory are 
fovind to compare favorably with the exact results of Floquet theory. 
Numerous numerical results are displayed as Brillouin diagrams for 
the first few Bragg orders. Moreover, explicit expressions for 
coupling coefficients, bandgap shifts and bandgap widths are derived 
for singly periodic media. Particular note is taken of phase speeding 
effects. 

The effects of multiharmonic periodicities on the control of 
feedback strength are investigated. It is found that with proper 
phasing the feedback strength becomes zero and the bandgap disappears. 
Coupling parameters are calculated for typical multiharmonic perio- 
dicities for the first three Bragg orders. 

For, odd Bragg orders, inverted bandgaps and phase slowing 
occur when the gain or loss of the media is modulated. Also average 
gain or loss affects the bandgap shape and the spatial or temporal 
growth or decay. Absolute instabilities are observed and expressions 
are derived for the instability frequencies, thresholds and growth 
rates. Under certain conditions, instabilities occur for structures 
with average loss. The results for the first and second Bragg orders 
are archetypical of all odd and even orders respectively. 


V 


Applications of the ECW theory to higher-order DFB filters 
involve such phenomena as transient propagation, effects of periodicity 
profiles and the relative coupling due to bovmdaries and periodicities. 
The calculation of higher-order DFB laser parameters shows that the 
mode spectrum is asymmetrically shifted and the threshold gain is 
greatly dependent upon the periodicity profile. Approximate threshold 
parameters are calculated for high and low gain and for all Bragg 
orders. In addition, application of the ECW theory to holographic 
gratings and beam propagation is made. 
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CHAPTER I 
INTRODUCTION 

This report investigates Bragg interactions in periodic media 
by using the example of electromagnetic waves propagating in spatially • 
periodic dielectrics. The main purposes of this report are to develop 
physically meaningful approximate methods for higher-order Bragg 
interactions, to show the mathematical foundations of the coupled 
waves theory and to give physically meaningful explanations for 
several previously unexplained mathematical results. Exact and 
approximate theories are compared numerically and applications are 
made to both bounded and unbounded media, to lossless and lossy media 
and to both passive and active media. Examples are given which 
correspond both to wave -packets in space and time and to the steady- 
state response of plane waves. 

The history of wave propagation in periodic media started with 

1 2 
Mathieu's equation in 1868 and subsequent generalizations by Floquet 

3 

and Hill in the 1880' s. Although Mathieu's equation had its origin in 
problems associated with elliptical boundaries, we will also show its 
connection to wave propagation in periodic media. This latter problem 

4 

was first considered by Lord Rayleigh in 1887 . He considered the 

effect of periodic density variations upon the propagation of waves on 

a string. In the early 1900' s a different, more physical, approach 

was taken by Sir William Bragg. He derived the necessary spatial 

period for constructive reflection of X-rays by crystals. These ideas 

44 

were formalized for quantum mechanical applications by Bloch in 

5 

1928. Two books in the 1940' s, one by McLachlan and the other by 
Brillouin,^ summarized previous work with Mathieu fvinctions and 


2 


with waves in periodic media. The books also provide useful 
bibliographies. 

While most of the above work was concerned with exact 
solutions of differential equations, a second, independent approach 
was taken in the late 1940' s and early 1950' s. This approach 
stressed the physical concept of wave coupling by periodic pertur- 
bations. A mechanical device demonstrated this effect in 1949 by 

coupling torsional energy between two bicycle wheels which were 

7 9 

periodically loaded with magnets. In 1953 Pierce used energy 
considerations to formulate what is now known as the coupled waves 
approach or the coupled mode theory. This approach has been 
popular because of its simplicity and intuitive appeal. Summaries 
of the coupled waves approach are given in texts by Pierce^^ and 
others. 18,26 last twenty years the coupled waves 

approach has been successfully used in such diverse areas as holo- 

gram diffraction, >•»>•>•> waveguide coupling, 

12,14,16,19,20 ^.17,74 

' * ' ’ traveling -wave tubes, ' parametric 

devices, X-ray diffraction,^^ distributed feedback (DFB) 
Users. ' 4 - 30 . 35, 36. 58, 70 others. 

37 

Extensive bibliographies on recent applications in optics and elec- 
tromagnetics are given in the references. We note that 
the telegrapher's equations which were developed before the coupled 
waves theory are of the coupled waves form. These equations are 
not approximations, however, since they exactly describe one- 
dimensional transmission line problems. 




The Floquet theory, which originated in the study of ordinary 
differential equations with periodic coefficients, has also been useful 
in the study of electromagnetic waves in periodic media. Although 
this theory is more cumbersome than the coupled waves approach, 
it provides an exact numerical solution. Extensions of the theory 


to include partial differential equations and finite length media have 
44 45 

also been made. ' Applications of Floquet theory to electro- 


magnetic waves have been made in the areas of traveling -wave 

anteimas,^® space-time periodic media, integrated optics, 

32 , 58 , 63 , 64 8,60 , 11 , 22 , 40 , 57 , 61 , 

' ' ' corrugated structures, ' and others. » » » • • 


62 , 75 


Other exact methods that are used in plane- stratified material, 
such as the matrix method^^ or the method of invariant imbedding,® 
will not be used here. 

The second chapter of this report contains the derivation of 
the first Bragg order coupled waves equations and the Floquet solution 
for electromagnetic waves in longitudinally periodic media. A simple 
explanation for phase speeding is given and the connection between 
the Floquet and coupled waves theory is explained. The dispersion 
relation is fotind as well as all pertinent coupling parameters. 

Several Brillouin diagrams illustrate physical principles and com- 
pare the approximate and exact theories. 

The primary purpose of the third chapter is to extend the 
coupled waves concept to higher Bragg orders. The resulting ex- 
tended coupled waves (ECW) equations provide explicit dispersion 
and coupling information for every Bragg interaction. Numerical 



4 




examples again illustrate the results of both approxiniate and exact 
theories. A section is devoted to effects caused by perturbations 
of several frequencies and the resulting disappearance of bandgaps. 

Periodic media with loss or gain is covered in chapter four. 
Inverting and non-inverting bandgaps are fovmd which depend upon 
Bragg order and coupling type. Index and gain/loss coupling are 
both considered. The effect of the periodicity upon average gain or 
loss near Bragg resonance is noted. 

The fifth chapter discusses the stability of active periodic 
media and gives explicit values for instability frequencies and 
thresholds at all Bragg orders. The stability criterion also speci- 
fies the correct root of the dispersion relation. 

Several applications of the ECW theory are given in successive 
sections of chapter six. The reflection and transmission of transients 
are discussed and demonstrated with detailed numerical examples. 

The extension of previous work to higher order DFB lasers is 
briefly covered. Diffraction efficiencies can be found when the ECW 
theory is applied to holographic gratings. Finally, the case of 
beam propagation in longitudinally and transversely periodic media 
is outlined. 

Conclusions of this report are given in chapter seven. 
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NX = 2 A sin9 (N=l,2,...) (2A. 1 ) 

where X = Zir/k = wavelength of plane wave 

A = 2ir/K = spatial period of structure 
N = Bragg order 

and where the velocity is assxxmed to be that of free space. The cases 
N:>2are referred to as higher-order Bragg interactions. For media 
with relative dielectric constant e, we restate the result as 

i 

ke^/K = N/2 (N=l,2,... ) (2A.2) 


for normally incident waves. Note that Bragg's Law does not account 
for the reflected wave amplitude, for the type of periodicity present 
or for the effect of slight variations of the wavenumber k from the 
value given by Bragg's Law. 

The latter effect is called phase mismatch and can be con- 
sidered in a semi-quantitative way by the use of Fig. 2. 2. This figure 
shows an incident wave I which is reflected from a three layer 
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b) 


c) 
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, 2 The effect of phase mismatch upon the reflected 
. Zero phase mismatch a) indicates 
2 

k = N/2. 


wave 
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periodic medium to form a total reflected wave R. The wave R is a 
phasor sum of four sub- reflections each of which has phases relative 
to the first sub- reflection at the incident phase. Multiple reflections 
are ignored in this simple model. Parts a)-d) of Fig. 2.2 show the 
relation of the strength of the reflected wave R to the phase mismatch. 
Fig. 2. 2a shows the constructive interference at the exact Bragg 
condition (i. e. zero phase mismatch) that produces large R. It is 
apparent that there is a considerable reflected wave R for slight 
phase mismatch (Fig. 2. 2b) whereas large phase mismatch will pro- 
duce small R (Fig. 2. 2d). 

Therefore, from simple wave interference arguments, one 
can deduce Bragg's Law and the qualitative effects of phase mismatch. 
Other theories are needed, however, to accoxmt for wave amplitudes 
and the effect of the form of the periodicity. 

I 

B. Coupled Waves Approach 

1. TEM Waves in Passive Unbounded Media 

Consider the case of a plane transverse electromagnetic (TEM) 

i 

wave that propagates in a longitudinally periodic unbounded medium 

as shown in Fig. 2.3. Assume a time variation of the form e 

S T 

starting with Maxwell's equations in a source-free, linear, iso- 
tropic region, we find in the frequency domain that 


7 X E (z, co) = 

itO H(z, to) 

(2.B1) 

7 X H (z, 05) = 

-itoe^e(z) E(z,to) 

(2.B2) 

7 • E (z, to) = 

0 

(2. B3) 

7 • H (z, to) = 

0 

(2. B4) 



-9- 
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Fig, 2; 3 TEM wave propagating in unbounded periodic media. 
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where CO = radian frequency 

free-space permeability 
free-space permittivity 

relative permittivity or relative dielectric 
constant 

and where E and H are the electric and magnetic field vectors. 
By combining (2. Bl), (2. B2) and (2. B3) we find the wave equation 
for no transverse variation 


to = 

^o = 

e = 
o 

e(z) = 


e(z) E(z) = 
-icot 


(2. B5) 

= transverse component of E(z, to) 


= to/c = 2ir/X = free-space wave number 
00 


dz 

where E(z, to) = E(z)e 

k 

e(z) = e[l + fp cos(pKz)J , f^ = 0, f j s i 

q s 1 is the perturbation. 

The periodic dielectric constant has been expanded in a Fourier cosine 
series. Assume that the electric field can be represented by just 
two waves, a forward wave F(z) and a backward wave 
^(z) which travel with positive and negative phase velocity along z. 
This assumption is intuitively appealing for q « 1 since these are 
the only two possible waves in the vinperturbed case. Thus, we con- 
sider the transverse electric field 


E(z) = F(z)e*^^ + B(z) e"^^^ 


(2. B6) 


where p is the longitudinal wavenumber. For first-order Bragg 
interactions, p/K 1/2. Then let 
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P - + A p = K/2 + Ap (2. B7) 

F(z) = F(z)e"^^^ B(z) = B(z)e"^^P"' (2. B8) 

Use equations (2. B5) - (2. B8) and the slowly varying approximation 


F(z) = F(z)e^^^^, B(z) = B(z)e 


F'f 

« ''F' 


( F" \ 

2 F 


B" 


J and 

(b'O 

« P 1 r. 

•^o V B / 

> 


primes denote differentiation with respect to z, to find 

C-p^F + 2ip F' + k^e F + ^ .5 ^ B]e^° 

•^ 0^0 2 

+[-p^ B - 2i p^ B' + e B + F]e‘^ 

n 2^ iSp I-, 2^ -,*i3P z 

= . 4n [e‘ J 

where the arguments in z have been dropped. The coupled waves 

±iPoZ 

equations are found by equating the coefficients of e to 


zero. 


15, 19 


The terms proportional to e 




(p = 2, 3, . 1 . ) are 


either ignored, termed non- synchronous or averaged over time and 
14 

considered zero. It will be shown later that these terms corre- 
spond to relatively unimportant coupling to other waves when p/K=‘l/2. 
The resulting coupled waves equations are. 



F'(z) - i 6 F(z) = i X B(z) 

(2. BIO) 


-B'(z) - i 6 B(z) = i X F{z) 
k^e-p^ 

(2.B11) 

where 

6 = ~~ 2 ~Q = phase mismatch 

(2.B12) 


^ ^ 4 ^ ^ ~ coupling coefficient 

(2.B13) 


- 12 - 


These equations agree in form with those derived elsewhere. 

18 20 23 30 

* ' ’ The equations account for both the wave amplitudes 

and the phase mismatch as well as the interaction of the waves with 
the fvindamental Fourier component of the dielectric periodicity. 

When 6 = 0, the waves F(z) and B(z) are coupled only through the 
perturbation T] while the change in amplitude of one wave is propor- 
tional to the amplitude of the other wave. For zero perturbation, 
the wavenumber of F(z) and B(z) becomes equal to the phase mismatch. 


2. Coupled Waves Dispersion Relation 

By differentiating the coupled waves equations, a wave equation 


[- 


dz'^ 


+ (6^-x^) 


][£ia ■ 


(2.B14) 


is constructed. Assuming a solution of the form the dispersion 

relation is found to be 


Ap = 



(2. B15) 


The approximations for 6 and y are 

1 1 
6 k e^-B = Ak 
^o 

X =“ xlk=k = ’I 
o 

'1. 

when Ap « K and where k^ e^/K = j and p^ = K/2. 
This produces the dispersion relation 



(2. B16) 
(2.B17) 


(2. B18) 


The following properties are evident for real e and r|. 
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1 . 

1. Waves propagate without decay for Ake VKl>|ri/8 | 
which correspond to passband regions. 

2. Stopbands or bandgaps occur for imaginary Ap/K or for 

1 . 

I Ake^ /K I <|ri/8 I . Here waves decay. 

3. The coupling of waves is maximum at the bandgap 
center (Ak e^/K = 0) where Ap/K = iri/8. 

The dispersion relation (2. B18) is plotted as a Brillouin 
diagram in Fig. 2.4 which clearly shows the regions of interest, 
namely the stopbands (ellipse) and passbands (part of hyperbola). 

Note that the coupled waves analysis gives much more in- 
formation than Bragg's Law. However, the coupled waves approach 

is only valid arovind the first Bragg interaction and does not describe 

41 42 

important wave interactions at higher Bragg orders. ' Further- 
more, we have assumed a solution that is based on only two waves, 

F{z) and B(z). This is strictly valid only as T] -♦ 0 when F(z) and 
B(z) are the two eigenmodes of the media. 

C. Floquet Solution 

1. TEM Waves in Passive Unbounded Media 
An exact solution to the wave equation (2. B5) may be con- 
structed through the application of Floquet' s theorem. A form of 
this theorem is stated as follows: a linear differential equation with 
coefficients periodic in z with a period A has a solution E(z) with the prop- 
erty E(z+A.) = e^^^E(z) where pis the fundamental wavenumber and A is 

x(3 z 

the fundamental period; define ?5(z) such that E = e ^ 0(z); then 

<<(z) is periodic in z since d(z+A) = d(z). A proof of the theorem 




-IS- 


IS given in 


several references.^* 


We now expand the fxmction ^(z), mentioned above, in a 
3 43-47 

Fourier series. * The resultant expansion for E(z) is 


E(z) = E a e 


i(p+nK)z 


n=-oo 


(2. Cl) 


which is the constructed solution for the wave equation. This solution 
is made up of an infinite number of space harmonics, a^, of order n, 
which propagate with longitudinal wave numbers (p+nK) (see Fig. 2.6). 
Substitute the solution (2. Cl) into the wave equation (2. B5) to find 


- (p+nK) ] a + 


E f 


n 2 p=-oo P 


I gi(P+nK)z ^ Q 


(2.C2) 


(n = 0, ±1, ±2,...) 


where we have defined f^ = f Algebraic manipulations transform 
this result to 


Da + Ea f = 0 

n n p =-«5 n-p p 


where D s -fl - 1 

n Ti L k^e J 


(2. C3) 


(2. C4) 


In matrix form this can be rewritten as 


D = 0 


where 


(2.C5) 


a 


(2. C6) 
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D = 


°-2 ^2 ^3 ^4 




^2 ^1 ^0 ^1 ^2 


‘^3 ■ ^2 ^1 ^1 ^1 


^4 ^3 ^2 ^1 °2 


The non-triviality condition for the matrix equation requires 


det p = 0 


which is the Floquet dispersion relation connecting p and k. For 

singly periodici media (i. e. = 0 for p s 2) an expression for the 

space harmonic ratios and for the dispersion relation can be found 

43 

in terms of a rapidly convergent continued fraction. 


2. Hill's Determinant 




suggested an alternative form for the dispersion 


relation which is eqmvalent to (2. C8). The derivation for the case 

42 

under consideration has been given in a previous report with the 
following resvilt. 

sin^(irp/K) = A(0) sin^( 7 r k /K) 


(2.C9) 


-17- 


where the elements of the Hill's determinant A(0) = det 1|‘^|| are given 


1 2 
- k e 

2 ^2 ,2 

p K - k e 


p = n 


p n 


2 

Figure 2. 5 shows the qualitative behavior of sin (ir |3/K) and the 
resvilting bandgaps where (3 is complex. 

3-5 

Several authors suggest an approximation for the infinite 
order determinant, 

1*+ i7? 9^777 

I 6 k€ 2 /K (Cl-4k e/K ] [2 -4k e/K ] 

[ 2 k^e T) f /K^]^ ) 

+ — s 5 ^ + . . . ? (2. C 

[3 -4k e/K ] ) 


wliich is valid for T]fp«l (p = 1, 2, 3, , . . 

Note the following points about the dispersion relation (2.C8 ) 
and its approximation (2.C11). 

1. This dispersion relation takes into account all Fourier 
components f^ of the perturbation and is valid for all 
k and p. 

1 . 

2. Away from the Bragg interaction regions ke^/K ^ N/2 

1 

(N = 1,2,...), A(O)'^ 1 as T) -• 0. Hence k ± p. Thus, 
the periodic medium has little effect upon wave propaga- 
tion in the passbands as T] -* 0. This was expected from 
Bragg's Law and the coupled waves approach. 






3. When sin^ (irp/K) > 1 or sin^(iT/3/K) <0, |3 is complex. 

From (2. C9) and Fig\ire 2.5, this occurs for k e^/K^N/2 
(N = 1,2,3,...). Hence, the dispersion relation agrees 
qualitatively with the position of the bandgaps predicted 
by both Bragg's Law and the coupled waves approach. 

Note that the Floquet solution provides all of the information 
of the coupled waves equations with the additional advantage of being 
an exact solution that is valid for all (3 and k. However, the Floquet 
solution lacks the intuitive appeal and simplicity of the coupled 
waves approach. In particular, one has to use an approximation to 
find the approximate bandgap placement. In addition, a tr\mcation 
of an infinite determinant must be performed. 

D. Relation of Coupled Waves Solution and Floquet Solution 

The recent interest^ 51, 52 relation of the coupled 

waves solution and the Floquet solution is important for two reasons. 

First, it gives a more rigorous mathematical foundation for the 

coupled waves theory. Second, it can be used to expand the coupled 

53 

waves approach to higher -order Bragg interactions. This will be 
shown in Chapter III. 

Consider equation (2. C3) for the case of cosinusoidal per- 
turbations. 


D a 
n n 


+ a + 
n+l 


n-1 


0 


D 

n 


ifi -ie±?K)i~| 
I L k^c J 


(n = 0. ±1, t 2, ...) 


(2.D1) 

(2.D2) 
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For T) « 1 


k e^/K -» ± (p/K+n) (n = 0, ±1, ±2, . . . ) 


The resulting Brillouin diagram is shown in Fig. 2. 6 to consist of an 

infinite number of space harmonics a^ which propagate as 

From Bragg's Law, it is known that important interactions 

1 

between F(z) and B(z) occur at k e^/K P/K j where the ± sign 
comes from considering waves of both positive and negative phase 
velocities. Therefore, consider only the a ^ , a^, a” and a^^ space 
harmonics in the Floquet solution for the electric field 


E(z) = a j 


-i(p-K)z + iPz - -ipz . + i(P+K)z 

! ^ +a e +a e ^ +a,,e''^ 

o o +1 


(2.D4) 


p - p^ + Ap = (±K/2 + A p) 


(2.D5) 


where the ± sign holds for space harmonics that intersect each other 
at p/K = ± i- 

- -iApz . + iApz, iPoZ 

E(z) = (a_j e ^ + a e )e 

^ ° / 
i F(z) 


, - -lApz + +iApz). -ip,^z 
+ (a e ^ +a,,e )e 

' o +1 ' 


i B(z) 


(2.D6) 


E(z) = F(z)e^Po2 ^ B(2)e-iPo2^ 


(2.D7) 




- 22 - 


This is exactly equivalent to the assumed form of the electric field 
in the coupled waves approach. The waves F(z) and B(z) are each 
the sum of two space harmonics which are in agreement with pre- ■ 
viously derived results. However, the terms (n=3, 4, 5, ... ) 

that were discarded in the coupled waves approach now have a clear 
meaning. These terms correspond to coupling to higher order 
space harmonics. Explicitly, the terms correspond to 

space harmonics a^^ ^^id a ^ which (along with a^j and j ) couple 

i 

to the adjacent Bragg intersections at k e^/K ^ +i, P/K ± 3/2, 

This truncated Floquet theory can also reproduce the coupled 
waves dispersion relation. Consider the first Bragg intersection 
at P/K Here only the a^ and a_j^ space harmonics are important. 

If we truncate the relation = 0 to include only these two space 
harmonics, we find the dispersion relation. 


where 



By evaluating the determinant we reproduce the results of the coupled 
waves theory, namely. 
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Similar results are found at p/K ^,-‘| if we use the a^ and a^^ space 
harmonics. 

We note the following. 

1. The coupled waves expressions for F{z) and B(z) 
are each the sum of two space harmonics. 

2. The coupled waves approach can be viewed as a trunca- 
tion of the Floquet solution. This shows that previously 
discarded terms correspond to coupling to higher-order 
space harmonics. 

3. The coupled waves approach shows that wave coupling 
in periodic media can be viewed as coupling between 
the intersecting pairs of waves that make up F(z) 
and B(z). 

E, Numerical Results 

The dispersion relations are numerically compared by using 
the Hill's determinant and the coupled waves approach. 

1. Limitations of Hill's Determinant 

The Hill's determinant dispersion relation 

2 2 - 

sin(irp/K) = A(0) sin (irke^/K) (2. El) 

is limited by the number of significant places used by the computer. 
In the case of the Univac 1108 this is 9 places for .complex calcula- 
tions. The smallest number that can be used which is larger than 
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- 8 

vinity is 1 + 10 . Thus, if 

sin^ (irP/K) ^ 1 + 10‘® (2. E2) 

at the first bragg order, then 

sin(irP^/K) cosh(trP^/K) + i cos(ttP^/K) sinh(TrP^/K) 

““ 1 + 10"® (2.E3) 

I 

where p = p^ + p^ . 

For -n-p^/K « 1 and ttP^/K = tt/ 2, the above ’relation is approx- 
inriated by 

1 + (iTpj)^/{2K^) 1 + 10"® (2.E4) 

through trigonometric expansions. The maximum value of p^^^/K is 

ti/ 8 from equation (2.D9). Thus, the Floquet numerical calculations 

-4 

should be limited to T) > 10 at the first Bragg order. In this report, 
_2 

T] ^ 10 to account for this and any other computer errors. Similar 
arguments limit the numerical calculations at higher Bragg orders. 


2. Brillouin Diagrams for Lossless Passive Media 

Each Brillouin diagram is a plot of normalized frequency 
1 1 


[k e^/K = U)e^/{cK)] versus normalized wavenumber [p/Kj . 

In many of the cases the diagrams are expanded around the Bragg 
interaction region where Ak and Ap replace k and p. 

Figure 2.7 illustrates the main features of Floquet theory for a co- 
sinusoidal perturbation (i. e. f =0forp?^ l)withr|= 1. The largest effect 

^ j. 

of the periodicity is in the vicinity of the Bragg wavenumber s ke ^/K = N/2 





Fig. 2,7 Brilloxiin diagram for first three Bragg orders when 
T| = 1, 0 using Floquet theory. Dotted lines show imaginary parts of 
P/K on separate scales. Dashed line represents the vinperturbed 
case. A cosinusoidal perturbation is assvuned. 




where the dispersion relation deviates from the unperturbed case 

(dashed line) and p becomes complex. Note that the bandgap is 

shifted towards the larger wavenumber and that this shift increases 

' 43 

with Bragg order. This causes phase speeding which is an in- 
crease in phase velocity (= 0)/P) due to an effective decrease in the 
dielectric constant. Physically, as the wave travels through the 
periodic medium, it speeds up and slows down with respect to the 

i 

unperturbed velocity c/e^. However, if we consider an average 
velocity (v) we find for singly periodic media 

1 ^ 

(v> = j J* ' ' • — dz (2, E5) 

0 A/e(l+*^ cos Kz) 

where A = 2 tt/K = spatial period. For T) < 1, expand the square 
root to find 

<v) = T — r (1-3- cosKz + cos^Kz + • • • )dz (2. E6) 

A €2 0 ^ ® 

<v> = [1+(3 it/16)ti^ + 0(Ti^)] ' (2.E7) 

2 

Hence, the phase speeding is accounted for by effects 0(T| ) for T| < 1. 

Figure 2. 8 is an expansion of the first-order Bragg interaction 
region of Figure 2. 7. The coupled waves dispersion relation (2.B19) 
is superimposed. Even for T| = 1, the coupled waves theory closely 
predicts the correct coupling coefficient as indicated by the maxi- 
mum value of P in the bandgap. However, the coupled waves theory 

N 

does not predict the bandgap shift or phase speeding for first-order 
Bragg interactions. 


K 
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0.3k \ 


/ 




. , . 0.05 r 


- 0.2 


-0.3 


T) = 1.0 

0.5 


Fig, 2,8 Brillouin diagram of first Bragg interaction with T] = 1. 

This compares Floquet theory (upper curve) with 
coupled waves theory (lower curve). Dotted and 
heavily dashed lines are imaginary parts of Ap/K, 
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The coupled waves result (lower curve of Figure 2.8) can be used 
for perturbations other than Tj = I by multiplying each scale by Ti. 

Figure 2, 9 shows similar results for r| = 1.0, 0. 1, 0.01 
from the Floquet theory. As T] decreases, the coupled waves theory 
becomes a better approximation to the Floquet theory. The curve 
for T] = 0. 01 is also the coupled waves curve for all three cases 
since the Floquet theory and the coupled waves theory are graphic- 
ally indistinguishable. 

Res\ilts of Figures 2. 8-2. 9 are summarized iri Table 2. 1. 



Coupled Wa’ 

s^es 

Floquet 


T) = l 

II 

o 

• 

>— • 

II 

o 

o 
1 — * 

ri = l 

T1 = 0. 1 

II 

o 

o 
) — • 

Coupling X/K 

0. 1250 

0.01250 

0. 001250 

0.1267 

0. 01250 

0.001249 

Bandgap Shift 
BGS 

o 

• 

o 

0. 0 

0. 0 

0. 04 

0. 0005 

<0.0001 

Bandgap Width 
W 

0. 250 

0.0250 

0. 00250 

0. 26 

0. 0255 

0. 00255 


Table 2. 1 Summary of coupled waves and Floquet theory at the first 

Bragg order. 


The-relative contributions of the different space harmonics 
are shown in Figure 2. 10. The upper and lower curves are the 
result of matrix truncation at 3x3 and 5x5 elements respectively. 
Each truncated matrix is centered around the matrix element A 

oo 

in equation (2. CIO). The differences are not large and the 5x5 matrix 
produces dispersion characteristics that are ^ 1^ different than those 
of the 19x19 matrices used in Figures 2. 7-2.9, 
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Fig. 2, 9 Brillomn diagram at first Bragg order for T] = 1. 0 (top 
curve), T| = 0, 1 (middle curve), and Ti = 0. 01 (bottom curve) using 
Floquet theory. Bottom curve also represents coupled waves 
solutions for all three cases. Note difference in scale for each 
case. Imaginary A(3/K values are the elliptical curves with 
Separate scale. 
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Fig, 2.10 The effect of limiting space harmonics in the Brillouin 

diagram using Floquet theory at the first Bragg interaction 
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F. Modifications and Comments 


1. Arbitrary Periodicities 


The coupled waves, equations were derived for 
00 

e(z) = eO+'H ^ ^ cos(pKz)) where only the f, term played a role 
p=0 P ^ 

inthe first Bragg order calculations. For generality, consider the 
expansion of a completely arbitrary (although smooth) e(z) as 


€(z) = e l+nS [f cos(pKz)+g sin(pKz)] 


= e ] 


p=0 ^ 


where f^ = 0 = g^ 


Tlf, « T,, 


flgi = 


The analysis is similar to the previous calculations which 
lead to the coupled waves equations (2. BlO-13). Following the 
identical procedure we find 

■ f'(z) - i 6F{z). = i B(z) (2 


where 


B'(z) - i 6B(z) = i x" F(z) 

1 2 o 2 

k e -p 

r o 


^ ■ 4 p 

The dispersion relation (2. D9) is modified by the substitutions 


n - (Tlf + r]^f 

x^ x" 


(2.F7) 


All previous Brillouin diagrams can be used with the above substi- 

— 4 - 

tutions for the periodicity given by (2. FI). The fact that x = (X ) » 

where the asterisk denotes complex conjugate, is a general, result 

1 8 

which holds for lossless systems. 


2. Corrugated Surfaces 

The previous results are strictly valid only for volume per- 
turbations. Similar results have been extended to surface perturba- 
tions or corrugations.^^’ The extension involves the assumption 
that the surface perturbation can be replaced by an equivalent volume 
perturbation. However, this assumption, known as the Rayleigh 

assumption, is valid only for Kd < 0.448 (where d = corrugation 

33 

depth and K = periodicity wavenumber) as shown by Millar. 
Physically this occurs because deep surface corrugations have 
proportionately less effect on the surface waves than do shallow 
corrugations. An exact Floquet analysis which solves the boundary- 


value problem has been given.- 


63, 64 


3, Comments 

This chapter has set the groundwork for the next calculations. 
In addition the differences and similarities of the Floquet and coupled 
waves theory were discussed. In particular the coupled waves theory 
was seen to be an approximation of the exact theory-, where certain 
space harmonics were retained and combined while others were 
discarded. A similar process will lead to descriptions at higher- 
order Bragg resonances in the succeeding chapters. 
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CHAPTER m 

HIGHER -ORDER BRAGG INTERACTIONS 

This chapter is an extension of the previous chapter to higher - 

order Bragg interactions (i. e. N ^ 2 in Bragg's Law). Present naan- 

made dielectric periodicities other than superlattices are limited to 

the order of A ~ 1000 A, Therefore, some applications in integrated 

optics reqxiire operation at higher -order Bragg interactions. Already 

40 7CL71 

higher order DFB lasers have been experimentally demonstrated. ' 
Other optical applications include couplers and filters. 

In section A, the coupled waves theory is extended to all Bragg 
orders for singly periodic media. Explicit expressions are given 
for all important parameters relating to the bandgap and coupling. 
Numerical examples are given in section B for the first three Bragg 
orders. Multiharmonic periodicities and a fourth-order numerical 
example are given in section C. An example of disappearing band- 
gaps is also shown. In some of the illustrative numerical examples, 
we will use large values of the perturbation that may not be physically 
realizable. However, the objective is to dramatize the effects of the 
perturbation on the Brillouin diagram. In addition, the extended 
coupled waves (ECW) examples are easily scaled for other values of r|. 

A.. Extended Coupled Waves (ECW) Theory 

1 . TEM Waves in Passive Unbounded Media 
By extending the assumptions made in Chapter II, we state 
the following assumptions for N^^ order Bragg interactions: 
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1. The most significant space harmonics are Fj(z) = a^e^^^ 

+ a,_^ ^-i(P-NK)z ^ e-iP%a_^ ,i(P+NK°z^ 

2. To provide cross -coupling between Fj(z) and Bj(z) we 

account also for the following pairs of space harmonics 

- + + - + . 

a_j^j ), (a^j^^2» ^^2 )'•••' ^N-1 ^ ''^hich are 

X 

slowly varying near k e^/K = N/2. In this way we accoiint 
for the intersecting space harmonics between Fj(z) and 
Bj(z) in the simplest possible manner. 

3. Self-coupling which occurs between Fj(z), Bj(z) and their 
adjacent space harmonics must also be included. 

4. Assume 11 « 1 although the theory may hold for T| 1 as 
in the first order case. 

5. All other space harmonics are ignored. 

The derivation is started by including the above mentioned 
space harmonics of the Floquet theory in the expression for the electric 
field 


0 


E(z) = E 


N+1 


n=-(N+l) 


n 


-i(p+nK)z ^ 2] a. 

_ ^ n 


(3.A1) 


n=0 


.th 


Near the N order Bragg interaction let 


p-p^ + Ap= (±NK/2+Ap) (3. A2) 

where the ± sign is used for space harmonics that intersect each other 
when p/K ^ 0, The electric field becomes 


/ 
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N*/2 


E(z) 


= S e^^*^^lexp[i(l -2n/N)p^z] 

n=-l *" 

+ ^2 ^ J. a - -lAPz Jexp - [i(l -2n/N)p.z] 

n=-l LN-n 

I 


where 


N*/2 = 


N - 2 


N - 1 


for N even 


for N odd 


Co)- 


1 for N = 1 


for N even 


0 for N odd 


This can be rewritten as 


N*/2 


where 


E(z) ' = S -[Fj_ 2 ^^j^(z) exp[i(l-2n/N)p^z] 
n =-l 

+ ®l-2n/N^.^^ exp[-i{l-2n/N)P^z]} + S(z) 


F (z) - a" + a'*' e^^^^ 

^l-2n/N' ^ ~ n-N ® ^ ^n ® 


, - . + lApz , 

L _ = a.. e + a e 

l-2n/N ' N-n -n 


c/ \ - -iApz + iApz rr / s T 3 / \ 

S(z) = ® + ^N/2 ® = ^o^^^ = ®o^^^ 


Substitute the assumed form of the solution (3. A4) into the wave 
equation (2, B5) and use the slowly varying approximation as before 
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We drop the arguments in z for simplicity. 

{ - (l-2„/N)"p2Fj_2/^/j,t2i(l-2n/N)p^F'_2^/j^}e'<*-2"/N)Po- 

+ {- (l-2n/N)^„2 . 2j,i.2„/„,p^B,_^^^J^-i(l-2n/N)p„. 


+ k^e I: 


-2„/N ^ B. e-^<> p„z^ 


^¥^(b.- 2„/K .-ia-^n/N,p„.| 


l-2n/N 
‘l-2n/N 


x{S fp[e^2pPo^/N^3-i2pPo^/NJj^Q 


(n = -1. 0, 1, 2 ) 


(3.A5) 


which is analogous to (2. B9). Next, for simplicity we limit the 
Fourier coefficients of the periodicity such that ^ for p ?? 1 
(i. e» singly periodic media). Equate the coefficients of 
to zero for each n to get the following N+3 coupled equations. 


(k^c - P^^Fj + 2i p^ F,' 


i2 

_ - k er] 


^^l-2Ag'^^l+2yN^ 


[k^e-(l-2/K)2p^^]Fj.2^+2i(l-2/M)pF,'.2^ = -k^,Fj+F, 


2 

(J)k"ss = -(J)-^<r-2/N+B2/N) > (3.^) 




_ -k^et) 


(B 


(k^€-p^)Bj + 2iP„Bj' 

[k^e-(l+2A^fp^]Bj^2A^-2i(l+2/l^)P„Bj^2^ = B, 


1 -2/N'^® 1 +2/N 
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i 

We define interaction as meaning the region in the (ke^, P) plane near 
the intersection of two space harmonics. We again note that higher - 
order interactions refer to interactions with frequency greater than 
that of the fundamental Bragg frequency defined by ke^/K = j. Thus 
each wave, F(z)^ and B(z)j^ 2 defines an interaction. The 

coupling diagram. Figure 3. 1, and the above set of coupled equations 
show that Fj(z) and Bj(z) are coupled through N intervening inter- 
actions. Since each coupling is proportional to r|, the cross -coupling 

between Fj(z) and Bj^(z) should be proportional to T]^. The adjacent 

2 

interactions to Fj(z) and Bj(z) contribute terms of order T] to the 
phase mismatch and hence both the F(z)j^ 2 /jvj 

need to be retained to obtain correct self-coupling. The F(z)j^ 2 /jvj 
and terms will not contribute to the cross -coupling except 

to order This small contribution is ignored. 

To show the above statements mathematically, apply the 
slowly varying assumption and solve the inner N-1 and the outer 

two equations of (3, A6) to get all waves ^(^)i_2n/N’ ^^^^l-2n/N 

(n = -1, +1, 2, 3, . . . ) in terms of Fj (z) and B ^ (z). The outer two equa- 
tions can be solved trivially for ^^^^i+2/N* 

We begin by using matrix manipulations on the inner N-1 
equations of (3.A6) 

Fi- 2/N = det||A||/det|lcll j 

®1-2/N 


= det l|G|!/det ||cl 


(3.A7) 



where 


C C- 


-39- 


C C. 


(3.A8) 


^2 


^2 C 

C C 


(3.A9) 


C 0 


G = 


(3.A10) 


^2 ^ 


C = € - (l-2n/N)^ 

n ' o 


C =Tike/2 


= Aj = -k^eriFj^ 
Cj = A^ = -k^eriBj^ 


Instead of solving (3. A7) exactly, we approximate the determinants for 


T] « 1 when N ^ 2. Since C « C , we find 
' n’ 

i 1 1 1 1 "T ■ f * ^2 . * — 1“', I ^ t 


det lldl f'"(n) - (^2n/N)^ ] 


1 
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(N-l)/2 


jy f (n) for N odd 
n=l 


where 


■* f^(n) 


N/E 

TT f(n) for N even 


n=l 

v< 1 


for N = 1 


Approximate det A and det jc in a similar manner. 


detllAll ^ TT*f^(n) +(-l)^A»C^"^ (3.A12) 


det||G|l IT* f^(n) + (-1)^ G, C 


(3.A13) 


The solution for F(z)j^_ 2 yj^ becomes 


vN 


*■'==>1-2^ = -C^ + 


1 . (-1) AjC 


N-2 


TT*f^n) 


(3.A14) 


.th 


Near the N Bragg interaction we approximate the free -space wave- 
number by 


k^€ - e = = (NK/2)^ . 


The expression for F(z)j 2/N approximated by 


-(Ti/2)F (z) {-k^e 1^/2^-^ B (z) 


l-(l-2/N)‘' TT f (n) 

Extend the above formulation to include the N=1 case by 
introducing the symbol 


(2.A15) 


} 
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_ TiNFj(z) (n/2)^"^ Bj(z) 

^ ^1“2/N b( 1-1/N) I I * {4n(n-N)/N^}^ 

r 0 N = 1 


where Cjj = 


N ^ 2 


Similarly for B(z)j^ 2/N the analogous result. 


B(z) 


1-2/N 


-Cj^T)NBj{z) (-D^'^Mri/2)^'^ Fj(z) 

8(1 -1/N) ]T *{4n(n-N)/N^}^ 


The expressions for F(z)j^^ 2 ^j^ B(z)j^ 2 ^'j^ trivially fotxnd from 

(3.A6) under the stated approximations as 


F(z) 


B(z) 


1+2/N 


1+2/N 


-Cj^(ti/2)Fj(z) Cj^nNFj(z) 


1 -(1+2/N)' 


8(1 + 1/N) 


-Cj^(n/2)Bj(z) Cj^riNBj(z) 


1- (1+2/N)' 


8(1+1/N) 


Expressions (3, A16-3, A19) are substituted into the following two 
equations from (3.A6) for F^(z) and B^(z) 




(k^e-p^^B,(z)+2ip^Bj'(z) = ;£^(B(z)j_2/n+B(^>1+2/n' 


Upon rearrangement the following extended coupled wave (ECW) 
equations are formed, 

F/(z)-i6^Fj(z) = iXj^Bj(z) ( 

-Bj'(z)-i Bj(z) = iXj^ Fj(z) 


(3. A23) 





(3.A24) 


(3.A25) 


N = Bragg order 

= nK/2 

The ECW equations agree in form with other work in diffraction 
22 68 69 

or holographic gratings ' ' where the periodicity is perpendiciilar 

to the propagation direction. However, 'this is the first time that 
analytic expressions have been derived for propagation in longitudinally 
varying media where the coupling, bandgap shift and bandgap width 
are given explicitly for all Bragg orders. 

Table 3. 1 presents numerical factors found in the ECW equa- 
tions for the first five Bragg orders. 


N 

]]^{4n(n-N)/N^}^ 

g(N) 

= (4/N)(l-l/N) 

1-Cj^(n/2)^{N^/[2(N^-1)]} 

1 

1 


0 

1 

2 

1 


1 

1 - I (n/2)^ 

3 

(8/9)^ 


8/9 

I 

4 

(3/4)^ 


3/4 


5 

(16/25)^ 


16/2^ 



Table 3. 1 Numerical values for factors in ECW equations. 
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2. ECW Dispersion Relations 

The dispersion relation is derived from the ECW equations 

. 6 

as in the first order case. 

As.. /76:,v2 ;y\, .2 


^ -Ji^f-C^f 


(3.A26) 


The following approximations may be made in the interaction region 


when A 6^^ « K. 




1 , c 2 .^2 

Ak.. 

^ 16S^(N^-1) 




2 

: i 

(>i 

,N+1 , .1 


Xn 


1 -ir- ko Ti- 


where 


^N+I 

1 J * [4n(n-N)/N^]^ 


i 

= 

k - p /e2 


P-Po 

Po = 

NK/2 


(3.A27) 

(3.A28) 


.th 


The N order dispersion relation is explicitly 


1 


N _ 


Ak.. e 2 


i 2 


N 


K 


K 


•).c 


Ak-.e2 2 ,..3 
N Ti N 


7 N*” 


N K 


16(N^-1) 1024(N^-1)^ 


2N 


.N+2 


N ) 

\'n’"{3n(n-N)/N^3^ / 


(3.A29) 


Note that for real e and r|, the maximum imaginary part of A p^^/K 


occurs when = 0. This defines the normalized coupling in terms 

of the bandgap behavior since. 
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lm{^} 

L K J max 


K 


.N 


N 


^N+2 -j-j-=:= [4n(n_isr)/N^}^ 


(3.A30) 


If Sj^/K = 0, then maximum coupling takes place at a wavenumber 

' jt. 

that is displaced from the exact Bragg condition ke^/K = N/2 for 


N ^ 2. This is referred to as bandgap shift (BGSj^) and is defined by 


BGSj^ = 


K 


5n=0 


32{N^-1) 


(3.A31) 


We note that this raises the Brillouin diagram and produces phase 
2 

speeding 0(ri ) as expected. The bandgap occurs whenever 
|6 j^/K| < 1xj^/K| which causes Ap^^/K to become imaginary. The 
bandgap width (W'j^) is defined by 


1_ r\^ __ N 

K 1=2N+1 {4n(n-N)/N2}2 


(3.A32) 


Table 3.2 summarizes some of the ECW results for the first 
five Bragg orders. < 


N 

Coupling 

X^/K 

Bandgap shift 
BGSj^ 

Bandgap width 
W 

1 

■ n 

0 

n 


8 


4 


2 

2 

2 

2 

. Jl_ 


n_ 


8 

12 

4 

-a 

243 Ti^ 

27 

243 p^ 


2048 

256 

1024 

! 

4 

2 


4 .. .. 

-3- 

2 r| 

2 T1^ 


9 

15 

^ '1 
9 

C 

1953125 Ti^ 

ne 2 

125 p 

1953125 p^ 

D 

18874368 

768 

9437184 


Table 3.2 Summary of the .main features of the ECW 
theory for the first five Bragg orders. 
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We note the following points about the ECW theory; 

1. The ECW derivation is an intmtively based theory 
that gives explicit values for x^/K, 

and for any Bragg order. 

2. The ECW theory predicts coupling coefficients that 

N 

are proportional to (T)/2) due to cross- coupling. 

This is expected because of the N interactions be- 
tween Fj(z) and Bj(z) where each interaction couples 
with strength T)/2 to adjacent interactions. The sign 
of the coupling coefficient alternates with Bragg 
order. 

3. The ECW theory predicts that maximum coupling 

occurs not at exact Bragg resonance but at a shifted 

frequency instead. This shift is proportional to 
2 

Nt| for N ^ 2 and is due to the self-coupling of 
Fj(z) and Bj(z). This bandgap shift accounts for 
the phase speeding effect that was first found in the 
Floquet results. 

N 

4. The bandgap width is proportional to {r\/2) . Hence 
for large N, only a small range of frequencies will 
cause significant coupling between Fj(z) and Bj(z). 
Also note that for large N, W^^ < 2 BGSj^^ so that the 
longitudinal wavenumber p may be real at exact Bragg 
resonance. Since bandgap width and coupling are 
proportional, it is impossible to attain large coupling 
and small bandgaps sipaultaneously. 
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j , 

5. A simple scaling rule exists such that if r) and e^/K 
are reduced by the same factor, then 6p^/K will be 
reduced by the identical factor. Thus ECW results 
for large perturbation can be directly applied to other 
perturbations. 

The extension to, sinusoidal perturbations is made in the same 

2 2 2 

manner as in the first order theory. Let '*’^2 ^ 

periodicities of the type rij cos(NKz) + sin(NKz). 

The relative magnitude of the waves F(z)j 2n/N' 

B(z)i 2n/N ~ 2,3,... N /2) can be found from equation (3. A6) in 
the same manner that F(z)j 2/N ^( 21)2 2 /N found. The 

results are in terms of Fj(z) and BjCz) v/hich are in turn related by 
the boundary conditions. We will not need these results now. 


B. Numerical Examples 

1. Second-Order Interaction 

The ECW approach uses the ten space harmonics shown in 
Fig. 3.2 for N = 2, The explicit dispersion relation is 


,2 


K 


K 


)-( 


5 - 

/ 6 ■ 576 ^ 


K 


(3.B1) 


and the phase mismatch and coupling are 

_ k^e (1 -t,^/6) - 

‘2 rr^ 

,2 2 

^2 8 p 


where 


p = K 


(3. B2) 


(3. B3) 



Fig. 3, 2 Coupling diagram which shows the ten relevant space 
harmonics that are used in the ECW theory for secon 
order interactions. 
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Note the change in sign of the coupling coefficient from that 
of the first-order theory. The coupling is provided by S(z) and the 
bandgap shift is accounted for by S(z) and F 2 (a). Phase 

speeding is the resvilt of the bandgap shift. 

Figure 3.3a,b, c shows a comparison of the Floquet and ECW 
theories at the second Bragg order for rj = 1.0, 0. 1 and 0.05 respec- 
tively. Although one would not expect the ECW theory to hold as r| -♦ 1, 
Figure 3.4 demonstrates the validity of the theory in this case. We 
note that for t| = 1, the Floquet theory predicts a slightly larger 
bandgap than does the ECW theory. However, in the practical case 
of T) ^ 0. 1, the two theories become graphically indistinguishable. 
Table 3.3 summarizes these results for second order Bragg inter- 
actions. 



ECW ' 

Floquet 


II 

• 

o 

T) = 0. 0 

T| = 0.05 

71=1 

• 

o 

II 

BS 

Coupling 

Xj/K 

. 125 

. 00125 

. 0003125 

. L33 

. 00125 

. 000311 

Bandgap Shift 
BGS2 

.0833 

.000833 

. 000208 

. 10 

.00080 

.00022 

Bandgap Width 

W 

^2 

. 250 

. 00250 

. 00625 

. 26 

. 0025 

. 00063 


Table 3.3 Summary of ECW and Floquet dispersion 
relations at the second Bragg order interaction. 

Figures 3. 4 and 3. 5 show the effect of truncating the ECW 
and Floquet theories. In Figure 3. 4, the outermost space harmonics 
^(aj , a^) and (a ^ , 3 . ^ ) of Figure 3. 2^ are not used for T) = 0. 1. In 
this case the coupling between Fj and Bj is given correctly since the 



























Fig. 3. 3 Brillouin diagram at second Bragg order for a) T) = 1* 0, b) T1 = 0, 1, 
and c) T) = 0. 05 showing Floquet and ECW resvilts. Imaginary Ap/K 
values are the elliptical curves with separate scale. Note the close 
agreement- of theories for f) ^ 0, 1, 
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Fig. 3. 4 The results of neglecting the effects of and B 2 
in the ECW theory at the second Bragg interaction. 


* 





A0 


K 

Fig. 3. 5 The effect of limiting space harmonics in the Brillouin 
diagram using Floquet theory at the second Bragg 
interaction. 




-52- 


cross-coupling is only dependent on S. However, the bandgap shift 

is not given correctly since S, and contribute terms of order 
2 

r| to the self-coupling and bandgap shift. Figure 3. 5 shows the effect 
of limiting the space Harmonics in the Floquet theory for T] = 1.0. 

The upper curve has a greatly enlarged bandgap width for a 3X3 matrix 
which includes the space harmonics a^ , a^ and a The lower curve 
adds the a^ and a_* space harmonics. The resulting 5x5 matrix 
produces resvilts that are graphically indistinguishable from the 19x19 
matrix used for all other Floquet results. 


2, Third-Order Interaction 

The ECW approach has twelve space harmonics ^(a a ^), 

(a'.a^ ), (a“,a 2 ’’), (a£,a^)) for N = 3. 
The results are similar to the N=2 case. The explicit dispersion 
relation is 




14499 T)^ 
1048576 


(3.B4) 


and the phase mismatch and coupling are 


5 


3 


X3 


k^e (1 - 9 t]^/64) - 

Tp 

^o 

81k^ e 
1024 p^ 


(3.B5) 
(3. B6) 


where P = 3K/2 

o 

Figure 3. 6a, b displays the results of ECW and Floquet calcu- 
lations for r| = 1. 0, 0. 5 at the third Bragg interaction. Note that for 
the first time the half width (=Wj^/2) is less than the bandgap shift 



- 0.3 - 0.2 - 0.1 0 0.1 0.2 0.3 


M 

K 



- 0 .© - 0 . 02 - 0. 01 0 0.01 0.02 0 .© 

M 

K 

Fig, 5, 6 Brillouin. diagram at third Bragg order for a) T) = i, 0 and 
b) T) = 0,5 showing Floquet and ECW results. Imaginary Ap/K values 
are the elliptical curves with separate scale. 
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(= BSGj^) which causes the entire bandgap to be above the exact Bragg 

wavenumber. This effect becomes more pronounced with increasing 

order and decreasing perturbation since ('n/2)^ whereas 

2 

BSGj^ ( t \/ 2 ) . Although the ECW theory does not approximate 
the Floquet theory for N = 3 as well as it does for N = 2 (compare 
Figures 3.3a and 3. 6a), the ECW approximation improves as r] 
decreases. This is an expected result since in the ECW theory all 
waves except Fj and B^ are assumed to be slowly varying. However, 
as N increases, more of the waves contribute to the coupling between 
Fj aind Bj. Table 3.4 summarizes the results of Figure 3.6. 



ECW 

Floquet 


T1 = 1 

ri = 0. 5 

T) = 1 

T) = 0. 5 

Coupling 

Xj/K 

.1186 

. 0149 

. 136 

.0155 

Bandgap Shift 
BGS^ 

. 105 

. 0264 

. 16 

.030 

Bandgap Width ' 

W 

^3 

.211 

. 0297 

.27 

. 033 


Table 3. 4 Summary of ECW and Floquet 
dispersion relations at the third Bragg 
order interaction. 

Figure 3.7 shows the effect of increasing the matrix order 
that is used in the Hill's determinant for the Floquet dispersion rela- 
tion (2. El). It is apparent that space harmonics of order n > 2N+1 
must be used to insure accuracy of the Floquet result. Graphs are 
shown for 7x7, 9x9 and 17x17 size matrices. As before, 19x19 size 

matrices have been used in the Floquet theory for other figures in this 
section. ' 









C. Multiharmonic Perturbations in ECW Theory 


Higher-order Bragg interactions are interactions in the region 

j. 

of k e^/K p/K ^ N/2 where N ^ 2. However, there can be competing 
.processes between different Fourier components, f^ cos (pKz), of 
multiharmonic (i. e. p = 1, 2, 3, ... ) periodicities for N ^ 2. In par- 
ticular we may find passbands where bandgaps existed in the singly 

periodic case. This latter fact has been of interest in solid-state 
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theory where the Saxon-Hunter theorem states that a forbidden level 
in an infinite lattice of pure type -A potentials and in a lattice of pure 
type-B potentials is also forbidden in any alloy containing both type A 
and type B potentials. As applied to our problem, this theorem im- 
plies that bandgaps formed by two dielectric periodicities f(Kz) and 
T (Kz) are just the bandgaps caused by T(Kz) and r'(Kz) separately. A 
number of counter-examples^^' have been given to the Saxon-Hunter 
theorem. The ECW approach will show explicitly the effect of multi- 
harmonic periodicities upon the bandgap. 

The extension of the ECW theory is straightforward. However, 

the results for arbitrary multiharmonic periodicities become cumber- 

1 

some. Several special cases will illustrate the general theory. 

From the previous sections we know that the waves between 
Fj(z) and Bj(z) are needed to couple energy from Fj(z) to Bj(z). We 
also know that the waves ^^^^1+2/N included 

to properly account for the bandgap shift. In addition, other waves 
that can couple significant energy between F^(z) and Bj(z) or that 
couple Fj(z) and Bj(z) to themselves have to be included. 
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Since the Fourier component [^T| cospK ] couples energy 
between intersections, p intersections apart with strength approxi- 
mately proportional to T] then coupling diagrams will again help 
to show the waves that should be included in the ECW theory. Figure 
3.8 demonstrates some of the possible couplings for N = 2. 

We list the possible couplings and strengths shown in Figure 


3. 8; 

2 

1. two first-order couplings (Tlfj) ; 

2. one second-order coupling “T] f^ ; 

3. one third- and one first-order coupling « (n f^(T| f^); 

4. one third-, one second- and one first-order coupling 


a (Tifj)(Ti f2)(Ti f^); 

5. one fourth- and one second-order coupling cc(r| f^)(T| f^, ) ; 

2 

6. one fourth- and two first-order couplings « ('nf^)(rifj) . 

The mathematical solution consists of writing a set of equa- 
tions analogous to (3. A6) where the right hand terms are augmented 
by all of the possible couplings by each Fourier component of the 
periodicity, f^. 

1. Order Bragg Interaction with fi and' /n 

Consider the periodicity made up of Fourier components 

[Tlfi cosKz + r|fj^cosNKz] at the Bragg order. We assume the 

non-trivial case where 0[(ri^,)^] ~ 0(r|fj^) so that contributions to the 

cross- coupling from the two Fourier components are of the same 

2 

order. The major bandgap shift is given by terms O [('nfj) ]. Figure 
3,9 is the coupling diagram which shows, the important space har- 
monies and couplings that are used in this case. 



Fig. 3,8 Coupling diagram that shows several of the possible 
cross -couplings from to and self -couplings 
for Fj when N = 2, The approximate magnitude of 
the coupling strangths are shown. 







-60- 


Since the coupling t) f only affects the equations for Fj(z) 

and Bj(z) and the coupling r) f ^ is used to couple Fj(z) to Bj(z) as 

‘well as the coupling Fj(z) Fj(z) and F ^ (z) F(z) ^^ 2 /^ 

“* F.(z), we can write by inspection from Figure 3.9 and equation 
* « 

(3.A6), 


-k^e 




1,2 
-k e 


(k‘e-p‘) F,+2ip^Fj= 


(J)k^cS = 


1 1 2 

1 . k e 


- (q ) 2 ^ ^1^^2/n'^®2/N^ 


) (3. Cl) 


-k^c 


(k"e- Po)Bl - 2iPoBi = '^[(^^i(Bi.2/i,+Bj^2/n)^'n^i] 


>2„2 


-k^G 


(k^-G - (l+2/N)X)Bi^2/N-2iO+2/N)(3^Bj^2/N = ^1^ J 


The only difference from the previous coupled equations (3. A6) is the 
change of notation q — Tlf^ and the addition of the T| and Tjfj^Fj 

terms. It is apparent from the N+3 coupled equations (3. Cl) that the 
bandgap shift will be due to rif j as before but that the coupling co- 
efficient will be the algebraic sum of terms involving Tif j and tlf^^. 
Solve the coupled equations as before to find 


I 


Fj(z) - i 6j^ 

Fl(z) = i Xj^ Bj(z) 

(3. C2) 

- Bj(z)-i 

Bi(z) = iXj^ Fj(z) 

(3. C3) 

k^{ 

^ 

l-C^(Tlfj/2)^C N^/2(N^-1)]}- 

(3.C4) 

- 
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Xn 

Ap 


1 2 
k e 


2 P. 


[ 


T]f 


/ i\N+l , ^ ,N 

N . 

2 33 


Tr"'{4n(n-N)/N^}^ 


jj (3. 


C5) 


N 


K 




(3. C6) 


Table 3.2 may be used to find the appropriate bandgap shift. The 


effective coupling coefficient is found by adding the values for Xj 


Xj^. Hence the results in equations (3. C2-3. C6) could be written 


down by inspection. 

An interesting case arises when Xj^ = 0* For example at the 
second Bragg order, X 2 = ^ implies 




K 


x,=o 


J. 

K 


Ak^ 


(nfi) 


K 


12 


(3.C7) 


which occurs when t] f ^ = (ri fj)‘^/2. Here the bandgap disappears but 

v2 


the phase speeding effect is still given by the (Tifj) /12 term in (3.C7). 
It is interesting to note that a periodic medium does not necessarily 
have bandgaps but it will still have the phase speeding property. 

In practical cases, if the effects of T)f ^ and r|fj^ cancel each other, 
other periodicities in the media may still cause bandgaps. 


Figure 3. 10a, b shows the ECW and Floquet results for 

v2 


± T)f 2 = ("nfj) /2 = 1/4 at the second Bragg order. In the first case 
the effects of T] f ^ and qf^ cancel, and the ECW theory predicts no 
bandgap. The Floquet theory shows a greatly diminished bandgap 
with a coupling coefficient ~ 00568. The Floquet theory 

bandgap can be reduced by an order of magnitude, and perhaps more 


2 2 

by slightly adjusting the ratio f^/tifj away from f^/r^fj" = 0. 5. For 



■q\2 ■ - 0-25 

/9n/K- 1.0 


/ ^ FLOQUET 


Fig, 3.10 Brillouin diagram at second Bragg order for a) r|f = ^0, 5, 
Tjf^ = 0,25 and b) rjfj = v 0. 5 , Tjf^ = -0,25 showing ECW and 
Floquet results for multiharmonic periodicities. Imaginary 
/K values are the elliptical curves with separate scale. 

Note the dependence of the bandgap upon the sign of T|f^, 
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2 

= 0.55, X 2 ^^ 0.0002. The second case shows the large 

bandgap that occurs when the effect of Tlfj and rif^ is additive. 

The vanishing of the bandgap is due to the relative phasing 

of f j and fj^ and is therefore dependent upon the symmetry of the 

periodicity. Analogous effects have been noted in the electronic 
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stopbands of crystals ' and in stability diagrams of transverse 

57 

magnetic (TM) wave propagation in periodic dielectrics. 


2. Fourth-Order Bragg Interaction with f^, f2 and f^ 

In the preceding case, the resulting coupling coefficient was 
a direct sum of the coupling coefficients of the two Fourier compo- 
nents. This is not always the case. To illustrate, consider a fourth- 

4 2 

order Bragg interaction where 0(r]fj) ~ 0(r|f2) ~ 0(11 f^). That is, 
the Fourier components rjfj, Ti f ^ and T| f ^ contribute values of the 
same order of magnitude (i. e. 0(11 f^) ) cross-coupling. The 

details of possible couplings and the relevant space harmonics are 
shown in the coupling diagram of Figure 3.11. Weaker couplings 

4 

such as those proportional to (T]fj) ri f ^ have been neglected. 

From Figure 3.11 and equation (3.A6) we write by inspection 
the following coupled equations. 





"t 










(k^e - 9/4 |3^)F3 + 3iS F 3 = 

o 2 o T 

(k^e - P|)F, + 2 ip„Fj ■ .= 

(k^e -ip|)Fi+ip^Fi 

k^e s = 

(k^e - i ( 3 ^)Bi-iP^b; 

o 2 “2 

(k^€ - p^)Bj - 2 ip^ b; = 

(k^e-9/4 p^)Bi-31p„Bi = 
O 2 ° 2 
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-_k_ej 2 _ .£ p . 

[f j(Fi+F| )+ S + Bj] 

^i^[fj(Fj+S) 4£2 Bi] 

,2 

[fj(Fi+Bi)+f2(Fj+Bj)] 

^^[fllBi+Sj+f^Fj 

[fjtBi+Bal+f^S+f^Fj] 


(f,Bj) 




) (3.C8) 


where the arguments with respect to z have been dropped. Solve the 
coupled equations as before to find 


F 3 

2 

= frifj 

^1 




(3. 

C9) 

Bl 

2 


B, 

2 

t 




(3. 

CIO) 

Fi 

2 


r^-iL f f 

L 16 12 

. inf.ni 
8 18 

f f^ + 

12 8 





■ 

^1 

f 

L 16 

ih 

7 r.3 

f^ _ H- f 
2 8 

/J 

(3. 

Cll) 

Bi 

2 


r ^-^f f 

L 16 ^T2 

.j]if f2 
8 12 

T| - 3 
"81. 

] 






[4h 

f +2L_ f 3_ 
2 8 1 

ni f f2. 

8 T2 


(3. 

C12) 

S 

= J^(F 
9 ' 


+ 2iii 
2 16 

f^.^f 
1 32 2 

■ \ 

(3. 

C13) 


These values are substituted into the equations for Fj(z) and Bj(z) 
in ( 3 . C8) to find the coupled equations for^Fj(z) and Bj(z). The 
resulting values for 6^ and are 




2 

o 




X4 


L 9 ^ 1 2 ^2 ^ 9 ^ 4 J 


( 3 . CM) 


( 3 . CIS) 


where = 2K and where smaller order terms are neglected. 

As expected, the above expressions reduce to the first-order 

case if ^ j = ^2 “ second-order case if = ^4 = 0, and the 

fourth-order case if ^2 ~ ^4 ~ However, the presence of the mixed 
2 

term f ^ means that the above residts cannot be obtained from 
simple superposition of the terms from each component. As in the 
previous case, a proper choice of f j, f^, and f^ would allow us to 
eliminate the bandgap. \ 

The different terms in the expression of X4 can be easily ex- 
plained as shown in Fig. 3 . 11 . There are four possible ways to 
achieve a fourth-order coupling with the three Fourier components 


f,, f_. and f .. T?iese are: 

F 2 ’ 4 

1. four couplings through f j ; 

2. any combination of one coupling through and two 
couplings through f^ ; 

3. one coupling through f^. 

The f^f-, term in equation ( 3 . CM) has a large numerical co- 
efficient because this coupling can occur in three different ways. 


- 0.06 - 0.04 - 0.02 0 0.02 0.04 0.06 0.08 

K 

Fig* 3, 12 Brillouin diagram at fourth Bragg order when = 0. 5, 

T)f 2 = 0,2, and f)f^ = 0,1 showing ECW and Floquet results. Imaginary 

9 

AB/K values are the elliptical curves with separate scale. 
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Fig. 3. 12 compares the ECW and the Floquet solution for 
the values Tlfj = 0. S, = 0.2 and rif^ = 0. 05. Note the close 
approximation of the ECW theory to the Floquet solution at the 
fourth Bragg order with several harmonic periodicities even when 
the perturbation is relatively large. 


D. Comments on ECW Theor 


The predictive abilities of the ECW theory have been s\xm- 
marized at the end of Section A. The numerical examples demonstrate 
the accuracy of the ECW theory at higher Bragg orders and for 
multiharmonic periodicities. The ECW theory is particularly attrac- 
tive for small perturbations where the Floquet theory is cumbersome 

and time consuming. This is relevant to present optical applications 

-3 

where often r) ^ 10 

The extension to higher Bragg orders opens up the possibility 
of using the many first-order resxilts of the literature at higher Bragg 
orders by the use of the proper and )^. Several extensions of 
this nature will be made in a later chapter. 

Two interesting results have been explained in this chapter. 
Firs^the phase speeding effect was observed to be a result of self- 
coupling and occurs regardless of the value of the higher frequency 
Fourier components f^. Second, the disappearing bandgaps were 
explained through simple interference effects of different periodicity 
frequencies. This opens up the possibility of controlling feedback 
through the relative phases of two harmonics since feedback or 
coupling strength is directly proportional to bandgap width. 


It is anticipated that a similar analysis could be applied to 
bovinded media, space-time periodic media and active media. The 
' latter case will be covered in detail in the next two chapters. 

53 

Parts of this chapter have been summarized elsewhere and 

the results have been confirmed through calculations that use the 

95 

method of multiple scales. 




r% fm r~^ 


-T^~'im-^i 
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CHAPTER IV 

COMPLEX PERIODIC MEDIA 

Although the study of wave propagation in periodic structures 

has been a popiilar subject since late last century, most of the 

analytical studies have dealt with simple (passive and lossless) 

23 3 0 58 

periodic media. The recent development of DEB lasers ’ ’ 

has extended the applicability of periodic structures to the case of 

complex (e = + ie^, T| = T)^ + iT)^) periodic media. In their analysis 

23 

of DEB lasers Kogelnik and Shank briefly discuss the Brillouin 
diagram near the first-order Bragg interaction when the gain co- 
efficient is modulated. Using Floquet theory, Wan^^' presents 

a more detailed analysis of the first-order Brillouin diagram for 
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active periodic media. Higher-order and multiply- re sonant 
DEB lasers have been mentioned very recently in the literature. 

In this chapter we investigate in detail both the ECW and the 
Floquet theory for complex media at the first few feragg orders. 
Analogous results will hold for higher even and odd orders. In section 
A we use the exact Floquet solution to plot the Brillouin diagram for 
several values of the perturbation and average dielectric constant. 

In section B we use the approximate but simple ECW solution to 
study in detail the changes in the Brillouin diagram for the first 
three Bragg orders in singly periodic media. A short discussion 
of m\iltiharmonic periodicities is given in section C, and section D 
summarizes the results of the chapter. 

As in the previous chapter, some of the parameters in the 
numerical examples will be larger than those that might be 


experimentally obtained. However, this exaggeration is again used 
to dramatize the periodicity effects. The ECW examples are easily 
scaled for other values of gain/loss and perturbation. 

A. Floquet Solution 

We consider the case of a TEM wave propagating along the 
z axis in an unbounded medium with a periodically modulated dielectric 
constant; 


e(z) = e + ie. + e E (r| 
r 1 r p 

+i'n.) f cospKz 
r Vip p 

where 



= Re{e} 

e. 

1 

= Im[e } 

Ti^ = Re{ri} 

'll 

= Im{r|3 

and Iq = 0, f ^ = 1, and 

br *pl ^ *■ 

P ^ 

1 - 


(4.A1) 


p = 1, 2, 3, . . . . Index coupling will refer to the case where = 0 
and gain/loss coupling will refer to the case where = 0. Note 

that gain (loss) media implies < 0 (e^ > 0) for e excitation. 

The Hill's determinant equation 


sin (irp/K) = A(0) sin (irke^/K) 


(4. A2) 


of chapter II still holds with minor changes in the matrix elements 
A : r , 

pn II p = n 


A 

pn 


-k^c 
p^K^-k^e 


T) f 


, e 

J?- Pl T 
2 e 


p n 


V. 

These elements make up the Hill's determinant A(0) = 


(4. A3) 


det ljA||. 
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Several examples of the Floquet solution for singly periodic 
media are given in the Brillouin diagram Figure 4. 1 for real fre- 
quencies. The example for simple periodic media (real e and T)) was 
found in chapter III and is given in curve 1 of Figure 4. la. There 
are two main characteristics. First, note the presence of bandgaps 
in frequency where ^ becomes complex. These are known as non- 
inverting bandgaps. Second, note the phase speeding that is a result 
of upward bandgap shift. 

The Brillouin diagram changes drastically for gain or loss 
coupling where the perturbation is imaginary. In curve 2 of Figure 
4, la we plotted the case for real frequency where e = and f] = ir^j. 
This corresponds to media with successive amplifying and lossy 
layers while the average gain (or loss) is zero. Note the inverting 
bandgaps (i. e. bandgaps in longitudinal wavenumber) at odd Bragg 
orders. The pattern of non-inverting and inverting bandgaps alternate 
with Bragg order. However, bandgaps of both types are shifted 
toward lower frequency. This latter effect can be understood from 
the average velocity which was given by 

<v) = + 0(n^] (4.A4) 

e2 

for jt) ) <1. It is apparent that phase speeding (positive bandgap 
shift) occurs for r) = and phase slowing (negative bandgap shift) 
occurs for T) = iq^. This also follows naturally from the definition 
of BGSj^ that was given in the previous chapter (3. A31). 

Let us now consider the case of complex media where only 
the real part of the dielectric constant is sinusoidally perturbed 



0.5 1.0 1.5 2.0 0 ±0.05 ±0. 10 ±0. 15 



Fig, 4, 1 Brillouin diagram for first three Bragg orders using Floquet 
theory and real frequency using: a) index coupling (curve 1) and gain/ 

loss coupling (curve 2) when = 0; b) similar results for finite 

e^/e^. The light dashed lines are the unperturbed values (Tl=0), 
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but average gain or loss is present (i. e. e = e^ + ie^, T) = r|^). The 


corresponding Brillomn diagram for real frequency is shown in 

'curve 1 of Figure 4. lb. For e. 4 0, we observe that Re{(3/K} is no 

longer a constant across the Bragg region and the effective spatial 

Z 5 76) 

gain or loss increases appreciably near the Bragg frequencies. ' 

The \inperturbed value (r|=0) is shown by the light dashed lines. 

When the gain or loss is sinusoidally perturbed, the spatial 

gain or loss diminishes near the odd-order Bragg resonances and is 

enhanced near the even-order resonances. This case is illustrated 

in curve 2 of Fig. 4. lb where there is an average imaginary dielectric 

constant and imaginary perturbation. Again the light dashed lines 

represent the unperturbed values. 

We see that the spatial gain can be either enhanced or 

diminished near Bragg resonances. This behavior depends upon 

Bragg order and perturbation or coupling type. 

1 . 

The sign of Im{p/K} and Im[k g^/K) (not shown) have not 
been specified because the correct root of the dispersion relation is 
dependent upon the stability of the wave. For passive media, > 0 
and the correct signs are chosen to indicate spatially and temporally 
decaying waves since no sources are present. Hence, the media 
does not support instabilities. However, instabilities may arise 
in active periodic media where < 0 for some (or all) z. In this 
case an analysis of wavepackets in the periodic medium must be 
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made. This case is taken up in the following chapter. For the 


present we will only examine monochromatic waves and will not 
specify the correct sign of the dispersion relation roots. 
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B, ECW Theory for Complex Media 
1. Analytic Expressions 

Since we again assume a solution to the wave equation of the 


form 


4 : 

E(z) = 

+ ®(==)l-2n/N exp[-i(l-2n/N)p^z ] + ) S(z) j (4.B1) 


where the symbols have the same meaning as in chapter III. This 
implies that the gain /loss and the perturbation of the media are 
small. Kence for a singly periodic effective dielectric constant 


e(z) = e^+i + e^{ri^+ iri^) cos Kz (4. B2) 

we assume 


\\\> hj . ^ 

However, as in the case of simple media, we might expect the ECW 
dispersion relations to be a good approximation to the Floquet results 
even as 

\\\> hj » I ^ • 


Following the identical analysis that led to the ECW expressions 

fVi 

(3. A20-3. A25) we find immediately that for the N Bragg order 

,2 






2P, 


(4.B3) 
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, ,»N+1 , 2 ^ /« . • « 
(-1) k (11^+111.) 

T T 1 


TT (4n(n-N)/N^J 


2i2 


(4. B4) 


where P^ = NK/2, Note that the coupling coefficient Xj^ is not 

affected by e^. The dispersion relation as before is 

^Pm\^ / / Xr’ 


m - m-m 


(4. B5) 


2, First-Order Interactions 


At the first order Bragg interaction we have the approxinriate 
expressions 

6 , 


K 


21 

K 


Ak, e 2 i c, 

1 r 1 

K ■^46 

r 

T| 4- iTl. 
r 


(4.B6) 


(4.B7) 


which give 



if - 

Ak. G 2 . e. V 2 

' *‘^1 S' ^ 

\ K 

/ 

K \4e / ^ ^ \ 

r 

. K / 



T|2 - 11.2 
'r 'i 

iri Ti. 
'r 'i 


, 

64 

32 


e. 

1 

2e 


(4.B8) 


where 


Akj = 


k - p /e 2 

^o r 


AP/ = P-P, 


Note that at the first Bragg order, only affects the phase mismatch 
6j,and T). only affects the coupling X| • index coupling, the 

maximum value of Im{APj/K} occurs exactly at Bragg resonance 
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1 

Ak, e 2 /K = 0. The value is 
1 r 



Hence, near Bragg resonance, the effect of spatial gain or loss and 

the perturbation add as the sum of their squares for index coupling. 

This enhancement of spatial gain or loss has been observed in the 

Floquet solution of Figure 4. lb and is shown in Figure 4.2a for the 

ECW solution for several values of e./c when Ti = 0. 1. Far from 

1 r ' 

the Bragg resonance Im{A(3j/K} = ± 6j/K, the unperturbed value. 

i 

The temporal gain or loss, Im{Akj e^/K], remains constant across 
the bandgap with value e^/4s^. 

2 

For gain or loss coupling, Xj changes sign and therefore 
A^j/K is real if €. = 0. This produces an inverted bandgap as shown 
in the previous Floquet results. The ECW results are shown in 
Figure 4.2b for several values of when T) = i 0. 1. Two classes 

of behavior appear that depend upon the sign of 

A = C|e./e^l - |t)./2|] . 

For A < 0, a bandgap of width 2 [t)^ /64 - /(4e^)]^ appears in 

Re{APj/K}, and Im{APj/K} = 0 at Bragg resonance. The two roots 
1 

of Im{Akj t^/K] are of opposite sign. For A > 0, there is no band- 

gap and the gain or loss coupling diminishes the average spatial 

gain or loss near Bragg resonance. This is due to the fact that the 
2 2 

coefficients of n- and e. have opposite sign in (4. B8). In this case, 

i, 

the two roots of Im{Akj e^/K} are of the same sign. 

Thus, the effect of finite e. can either enhance (t) - T)^) or 
reduce (r) = iflj) the effective spatial gain or loss near the first Bragg 
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Fig. 4,2 Brillouin diagram at first Bragg 
order using: a) index coupling; b) gain/loss 
coupling; and c), d) both couplings in the ECW 
theory. 


±0 005 
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resonance, For both cases of coupling, increased e./e tends to 

1 r 

mask the effect of the coupling or perturbation. 

Figures 4.2c, d show the effect of having complex coupling 

(T) = + i Ti^) in a media with no average gain or loss. As the ratio 

Ti./'n^ is increased, the bandgap region changes from the shape 

typical of index coupling to the shape typical of gain or loss coupling 

1. 

for Re{A(3 /k] and RefAk, e ^ /K}. The values of Im[A(3 /K} and 
1 1 r 1 

i 

Im{Akj £^/k} tend to peak in the Bragg interactiop region. 


3. Second-Order Interaction 

The second-order Bragg interaction and all even -order inter- 
actions produce very different dispersion relations from those of 
the first-order for gain or loss coupling. This was seen in the 
Floquet results of Figure 4, 1. The fchange is due to the term (i T)^) 
in the dispersion relation which takes on different signs depending 
upon the oddness or evenness of N, In the second-order Bragg region. 


the approximate relations are 

62 ■ Ak 2 G^^ (TI^ + iTl.)^ ie. p (Tl^+iTl/. 

K“K “ .12 ■'■2el:" 6 

r 


Xz 

■ ■ 

K - 


(Tl^ + iTl^)‘ 


(4.B9) 


(4.B10) 


where Ak_=k-(3 /e^ 
2 or 


P = K. 
o 


Note that the perturbation affects not only the coupling X 2 ' 

but also the phase mismatch 62 . In fact, T) now acts to modify the 

2 

gain or loss from e ^/(2 e^) to the effective value e^(l-T) / 6 )/( 2 e^). 

The Brilloxiin diagram for the' case of index coupling is shown 
in Figure 4. 3a for various values of Cj/Er 'with T) = 0.1. The result 
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is similar to the first-order case with the exception of the bandgap 

shift. Again, at the center of the bandgap, the effective gain or loss 

and the perturbation add as the sum of their squares. The spatial 

gain or loss is greatly enhanced near Bragg resonance when 

le./e^l « I'nl^ whereas the temporal gain or loss remains constant. 

In the case of gain or loss coupling, a difference occurs since 
2 ' 4 

X is proportional to Tl in (4. BIO) and the bandgap shift changes sign 
which causes the Brillouin diagram to be a mirror image of the index 
coupling case about the axis ~ That is, at the second Bragg 

order, the difference between index and gain or loss coupling is a 
matter only of positive or negative bandgap shift. This is demon- 
strated in Figure 4. 3b for r] = i 0. 1 and several values of 

2 

This dependence of the phase mismatch upon t) accounts for the 
phase speeding and slowing effects shown by previous Floquet re- 
sults. The temporal gain or loss is constant as in the index coupling 
case. 

1 

Figures 4. 3c,d are similar mirror images about Re[Al^e^/K} = 0 

for = 0 and various ratios 'n-Z'nj,* The symmetry is expected 

since Tl always enters as an even power in the dispersion relation. 

1 , 

From (4. B9-10) it is evident that Im{A p^/K] = 0 at Ak 2 = 0 

only when ] ri^ j = 1 ti^ ]. This corresponds to curve 5 and is similar 

in shape to the gain or loss coupling Brillouin diagram at the first 

1 

Bragg order. For Im{Ak 2 G^/K} = ± r|^T)j^/12, Im{AP 2 /K} = 0 for 
T)^ $ T)^. Thus for complex coupling, the perturbation will mask the 
spatial gain or loss for some frequency near Bragg resonance, while 
a large spatial gain or loss will occur for nearby frequencies. In 


X 
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Fig, 4. 3 Brillouiii diagram at second Bragg 
order using: a) index coupling; b) gain/loss 
coupling; and c), d) both couplings in the ECW 
theory. 


200 O' 
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all cases of both index and gain or loss coupling, the temporal gain 
or loss is peaked in the vicinity of Bragg resonance, 


4. Third -Order Interaction 

From the structure of the dispersion relation it is apparent 

that odd Bragg order interactions will exhibit similar Brilloxiin 

diagrams for index and gain or loss coupling. To exemplify this, 

consider the third Bragg order approximate parameters, 

1. 

6 Ak T) 2 on f i3e. 

■i “ - 256 (* B 


243 . . .3 

- 2048 


where Ak- ~ k - P /g 2 

3 or 

6 =3 K/2 


Figure 4, 4a, b demonstrates the dispersion characteristics 
for index and gain or loss coupling respectively. Except for the 
bandgap shift, the results are nearly identical to the first-order case. 
That is, the perturbation and gain or loss add at the center of the 
bandgap to enhance the spatial gain or loss for index coupling. In 
the case of gain or loss coupling the Brillouin diagram exhibits two 
distinct behaviors which depend upon strength of the gain or loss 
relative to the perturbation. For small average loss or gain, the 
effective spatial loss or gain is zero at the bandgap center. 

The case of both index and gain or loss coupling is shown in 


Figure 4,4c, d. In all third-order cases, the temporal gain or loss 
is peaked in the Bragg region just as for the first-order case. 


Fig, 4,4 Brilloiiin diagram at third Bragg 
order using: a) index coupling, b) gain/loss 
coupling and c/, d) both couplings in the ECW 
theory. 
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C, Multiharmonic Periodicities 

41 71 

Recent experimental work in higher-order DFB lasers ' 
c 70 

and in multiply -periodic DFB lasers make the results of multi- 
harmonic ECW theory timely. As in the passive -lossless case, ECW 
theory cannot give explicit dispersion relations for arbitrary Bragg 
order N, since the number of significant Fourier components (f^) 
have to be specified. However, the extension is straightforward, 
although laborious, for any specific case. 

Consider the case where fj and f^^ are the significant Fourier 
components of the periodicity and 


e(z) = + i + ®r ^^1^1 cos(Nkz +0)] 


(4. Cl) 


where rip = (t)^ + i Ti^)p , P =(1, N) and O(Tij^fj^). 

Following the derivation of the previous chapter we find the slightly 
modified ECW equations; 


where 




Fj'(z) -i6j^ Fj(z)=iXj^'^Bj(z) 


■Bj(z) - i Bj^(z) i 


]}-p 

t JN 1 - *- 2(N -I)-* 


2 

o 




2 jri fj.[cos9±isin9] (-l/^'^^(rijf 


(4.C2) 


(4.C3) 


± _ k e ; 'N N 
- 2(3, 


,N 


TT*{4n(n-N)/N'^3 


2i2 


(4. C4) 
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The resulting dispersion relation is 


Ap 


(^) = (1) 


Xn Xjsj 


K‘ 


(4. C5) 


By varying the phase' between the Fourier components, Xj^j can be 


varied while 6j^ remains constant. This allows control of the coupling 

1 


and bandgap width [= i(Xj^ ^ without change in the bandgap 

shift. 


For N = 2, the explicit dispersion relation is 


(^)‘ = [ 


A k_ e ^ 
2 r 


(Tlifi^ 


K 


12 


] ■ 16 0^2^2^ 


2 , ^’^l ^1^ 


-(Tljfj)^ll2f2COse] (4.C6) 


Figure 4, 5 shows two examples of multiharmonic Brillouin diagrams 

at the second Bragg order. In Figure 4. 5a, the index coupling case 

2 

is shown when (T|j^ f j ) /2 = “^ 2^2 “ ^ ~ bandgap width 

and I^^^P2^^^max vary as sin(0/2). Figure 4, 5b represents 

the case of index coupling with spatial wavenumber K and gain/loss 

coupling with wavenumber 2K. The result shows not only the 

Re{Ap 2 /K] typical of gain/loss coupling, but also the effective 

gain/loss typical of index coupling. In both cases shown, 

1 

Im{ Al<^e^2 /K} = 0 for RelAp^/K^* The results for lossy media are 
similar if « X 2 ^ "^ 2 ^^' 


D, Comments on Complex Periodic Media 

The extension of chapter III results to complex media has 
been mathematically straightforward for both the Floquet and the 
ECW results. However, inter esting.new features appeared in this 
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Fig, 4. 5 Brillouin diagram at second Bragg order using mviltiharmonic 
periodicities in the ECW theory: a) index coupling with relative phase 
0 between two periodicities; b) index and gain/loss coupling. Note 
Im {AkeJ /Kj = 0 for Re [Ap/Kj. 
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1. Even and odd Bragg orders produce non -inverting 
and inverting bandgaps, respectively, for gain/loss 
coupling. Index coupling always produces the usual 
non -inverting bandgaps. 

2. Phase speeding and positive bandgap shifts are pro- 
duced by index coupling while gain/loss coupling 
produces phase slowing and negative bandgap shifts, 

3. The average effective spatial gain or loss (i. e. 
Im[APj^/K}) is significantly enhanced or diminished 


near Bragg resonance whenever "H | . En- 


hancement or reduction depends upon coupling type 
and Bragg order. 

1 

4, Temporal gain /loss (i. e, Im{ Ak^ /K}) is either 

constant or is peaked at the bandgap center. 

5. .Multiharmonic periodicities offer great flexibility 
in changing the bandgap shape through phase and 


coupling variations. If the Fourier components f 


are dynamically generated, this opens the possibility 
of quickly controlling the feedback strength in DFB 
lasers or filters. Possible applications are in the 
areas of microwave and optical filters, switches 
and modulators. 

In this chapter, Floquet and ECW results were not directly 


compared because little or no graphical differences appear whenever 
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ItTjK \\\> Ni/®rl 

for the first few Bragg orders. The accuracy of the ECW resvilts is 
similar to that of the passive-lossless case. The only notable differ- 
ence occurred in the first order gain/loss coupling case where the 
ECW results' predict Im{A(3j/K] = 0 at exact Bragg resonance. The 
Floquet resxolt gives a small but finite value here for Im{Ap^/K]. 

Portions of this chapter have been summarized elsewhere. 
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CHAPTER V 

STABILITY OF BRAGG INTERACTIONS 
IN ACTIVE PERIODIC MEDIA 


The problem of wave stability was briefly mentioned in the 
previous chapter with regard to choosing the correct root of the 
dispersion relation for the frequency and longitudinal wavenumber. 

It is evident that the preceding monochromatic plane wave analysis 
is inadequate to predict the correct spatial and temporal behavior 
of waves in periodic media in all but the simplest (passive and loss- 
less) cases. Our purpose is to formulate the stability problem in 
such a way that the preceding ECW equations can lonambiguously 
define wave characteristics. 

The problem of wave stability was greatly clarified through 

72 

the original work of Sturrock “ which provided a method for distin- 
guishing between growing waves (unstable media) and decaying waves 
(stable media) through the dispersion relation, Sturrock further 
divided unstable media into convectively unstable and absolutely un- 

c 

stable media. The former referred to waves which grew and then 
decayed at a given point in space like waves in traveling -wave ampli- 
fiers. The latter referred to waves which grow everywhere in space 

as those waves found in backward wave oscillators. This work was 

73 

extended by Briggs and was formulated as a mathematical pre- 
scription which unambiguously classified waves according to their 
stability. A large number of electron-stream interactions demon- 
strated the usefulness of Briggs' jnethod in plasma physics. The 
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stability criterion was again extended to include time -varying media 

56 

by Cassedy who examined the stability' of parametric interactions 
due to space -time periodicities. 

In this chapter, the stability of time -independent active 
periodic media is studied at the first- and higher -order Bragg reso- 
nances. The stability classification depends only upon the dispersion 
relation or the Brillouin diagram and hence the ECW theory of the 
previous two chapters will be used. We will follow the procedure 
suggested by Briggs, This requires the formulation of the field 
response as a function of a source, localized in space and time. 

In section A the ECW equations are developed with sources 
present and the stability criterion is stated. The application of 
the criterion to the first two Bragg resonances is carried out in 
detail in section B, These resonances are the archetypes for all 
odd- and even-order Bragg interactions. Explicit values for insta- 
bility threshold, frequency and growth rate are given. Section C 
contains a brief explanation of the effects of complex coupling (i. e. 

T) = T] + iT].) and multiharmonic periodicities upon instability param- 
r 1 

eters. A short discussion of the results is given in section D, 

A, ECW Equations with Sources 

We carry out the usual manipulations with the Maxwell eqxia- 

tions that include a current density source J(z). The resulting in- 

6 7 

homogeneous wave equation for a TEM wave in longitudinally 
periodic media is 
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r + k^e(z)”|E(z, o)) = -iOJ/i J(z, uS} 
^ ° 


(5. Al) 


-00 

e(z) = € +ie. + e S(ti +iri.) f cos pKz 
' ' r X r n' 'r 'I'p p ^ 


p 1 P P 

Assume that E(z) is made up of the N+3 waves ^(^)j^2n/N' • • •> 

Bj(z), B(z)j^ 2 j^/j^ that are used in the ECW formulation. We then 
find the following N+3 ECW equations for singly periodic media at 
the Bragg order. 


2 




+ -i(l+2^)(3 z 


[k‘€-p„]F, + ap^Fj 


' - -Jsl£ll 






[k^e-p23Bj - 21P^b; = 


(5. A2) 

where the arguments with respect to z and t have been dropped and 
the primes denote differentiation with respect to z. The current 
densities J ^(z)j 2n/N portions of J(z) which are in phase 

synchronism with the rest of the equation. We accomplish this by 
trxmcation in wavenumber space such that 


ia+2/N)Pz 




Jl"^_2n/i^P:{:(l-2n/N)Po)=J(P+(l-2n/i^)Po)rect{[p:f{l-2n/N)p^]/Ap } (5. A3) 


(n = “1, 0, 1. 2, ... ) 
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where 


J(P) 


rect(P) 


Ap 


J* J(z)e dz 

Ti ipi<i 

\o lpl>i 

« p 


X 


As we shall see, the exact truncation details are unimportant since 
the current will only contribute to the field at wavenumber P = P^. 
Since we are particularly interested in the waves Fj(z) and Bj(z), 
assume that the current is slowly varying in wavenumber over the 
N+3 significant waves. Thus, we consider only the current contribu- 
tions Jj(z) and Jj(z) and ignore current contributions 0(ri). Using the 
above approximation and the other usual ECW approximations we find 
the coupled equations. 


Fj(z,ax)-i6j^Fj(z, w)=iXj^Bj(z,o)) J+(z, w)e'"P°"^ 


I tou 

-Bj(z, co) -i6j^Bj(z,w) =iXj^Fj(z,w) -^^.2 (z,w)e 


iPo^ 


(5.A4) 


These equations can be solved by taking Fourier transforms with 
respect to z and solving for Fj(Ap, w) and Bj(Ap, oo). The inverse 
transforms produce the time and space variation. 


Fj(z,t) 



n 

-00 


co[(Ap+6)jJ(Ap+p^,w)-xJi'(Ap-p^,to)}i^^^"‘^^^ 

D(Ap, w) 


dApdoJ 

(5. A5) 
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i/i 00 w[(6p+6)j“(Ap-f3 w)-xJwA6+6^.W)}e^^''"^^^ dApdoJ 

. i;?: 


■where , 2 ^2 . 

D(Ap,0J) =. (-^) - (|) + (| ) = 0 


(5. A6) 
(5.A7) 


is the dispersion relation and the radian frequency to is related to the 

i 

■wavenumber by to = c(k^+Ak)/e^ and S’S before. All sub- 

scripts N have been dropped for simplicity. Assume that the current 
is turned on at t = 0, hence J(z, t) = 0 for t < 0. 

It is apparent that the integrand of (5. A5-6), excluding the 
currents, is the Green’s ftanction for a periodic meditim in the ECW 

approximation. It is of similar form to an exact expression which 

75 

accounts for all space harmonics and current components. 

As pointed out by both Sturrock and Briggs, it is necessary 
to investigate wave packets in space and time due to localized sources. 
Therefore, the previous ECW equations (5, A4) with a source turned 
on at t = 0 are a useful start. The remaining problem is to find how 
Fj (z, t) and B^(z, t) behave asymptotically by choosing the contours 
of integration correctly and unambiguously. Note that the poles of 
the integrand, or the roots of D(Ap, co), determine the field response. 

For causality, it is necessary to carry out the frequency 
integration above all singularities. The integration contour is shown 
in Figure 5, la. Note that no roots appear in the region enclosed by 
the contour for t < 0, This contour is the rotated Bromwich contour 
that is used for Laplace transforms. The contour in the wavenumber 
plane is ambiguous since the roots of D(Ap, oj) may cross the real Ap 
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axis for different values of 0)(= OJ^+ico.). This difficulty has been 

*7 3 ' 

remedied by Briggs by use of the following physical criterion: 

for a source with sufficiently large temporal growth rate (i. e. large 
tOj = s > 0), the waves away from the source must spatially decay. 
Therefore, if the contour is chosen in the to plane to satisfy causality 
for iO. = s, the roots of D(Ap, o)) will determine the response for F(z, t) 
and B(z, t) and must represent spatially decaying waves. Figure 5, lb 
shows this condition. Since it is convenient to integrate along the 
axis U3^ = 0, we depress the contour around singularities down to the 
real coaxis (Figure 5. Ic), As cOj^ "* 0 for fixed co., the roots of 
D(Ap, co) will: 1) stay in the same half plane of A(3; 2) cross the real 
Ap axis; or 3) merge from opposite sides of the real Ap axis in pairs 
for some 00^ = a > 0. As shown by Briggs, the above cases lead re- 
spectively to: 1) decaying waves (stable media); 2) amplifying waves 
(convectively xxnstable media); or 3) time growing waves (absolutely 
unstable media). In the first case, the contour in the Ap plane along 
real Ap axis is not deformed (i, e. decaying roots) while in the second 
case the contour is deformed (i. e. amplifying roots) by root crossings 
(Figure 5, Id), In the third case, as Cl\ o for some CO^ = c»)^ (Figure 
5. le), the merging root behavior in the Ap plane (Figure 5. If) pre- 
vents the distortion of the contour since the integration must be 
carried out between the merging roots. The details of the integration 
are carried out in Appendix A for absolute instabilities with the 
proper contours shown in Figure A, 1. It can be shown (equation (A, 7) ) 
that the field varies as 




mm 


Roots of 
D(A /S,w) = 0 


mm 
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Decaying 
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Fig, 5,1 Integration contours in a) to plane for causality, b)toplane 
for t > 0 and time growing source, c) to plane with deformation to 
to^ axis, d) Ap plane v/hich corresponds to c) with changing roots, 
e)tO plane for absolute instability at tO^ = <y , f) Ap plane corre- 
sponding to e). 
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F(z.t) « 
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(5. A8) 


Note that the frequency and wavenumber of the instability are just 


, « I 

that of the merging roots (i. e. A|3 = A(3 , CO = W = CO^ + i a). Therefore, 


ot 


the waves grow in time as e . The absolute instability parameters 
are determined from the mapping of the dispersion equation roots in 


the Ap plane for variable to^ and different constant values of co^. 


B» Application of Stability Criteria to Bragg Resonances 

The Brillouin diagrams take on different characteristics for 

L . . _ 

even and odd Bragg orders in the case of index (T) = T)^) and gain 

(T| = i T)j) coupling. Thus, we will examine in detail the first and 
second Bragg resonances for singly periodic media since they are 
the archetypes of all even and odd Bragg order n interactions. 

The first-order Brillouin diagrams are repeated in Figure 
5,2a, b for the case of index and gain coupling respectively with 
various values of average gain or loss. Figures 5. 3 and 5.4 show 
a mapping of the roots of D(Ap, oj) in theAp plane as ca varies from 
large positive values to zero for several values of co . Normalized 

i. 

frequency Ake^/K is used instead of cti Figure 5, 3a shows the case 
of no average gain or loss (e. = 0) and index coupling which corre- 
sponds to curve 1 of Figure 5.2a, The merging root behavior is 

1 i 

noticed for Re{Ake^/K} = 11^/8 when Im{Ake^/K] = 0. Thus, the 
instability has no temporal growth and the medium is actually stable. 

i 

This is also noted from the fact that Im{Ake^/K3 = 0 for all Re[Ap^} 
and hence no instability is possible, ' Therefore, the physical notion 
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of the stability of passive periodic media is confirmed. Since there 

are no root crossings, the proper choice^ of the sign of Im{A(3/K} is 

that which indicates spatially decaying waves. Figure 5. 3b shows 

the case of index coupling with average loss (e> 0) which corresp&nds to 

curve 3 of Figure 5.2a. There are no root crossings or mergings \ 

X 

and hence the media is stable. The proper signs for Im{Ake^^/K} 

and Im[Ap/K} indicate temporally and spatially decaying outgoing 

waves. The case of index coupling with average gain (e^ < 0) which 

corresponds to curve 3 of Figure 5.2a is shown in Figure 5. 3c, In 

1 

all cases, similar curves exist for Re[Ake ^/K} < 0 which are mirror 

r 

images about the imaginary Ap/K axis. Note that the root crossing 
1 

occurs for Ake^/K = ie. /4e ± T| /8, Thus, absolute instabilities 

r 1 r 'r ' 

occur for index coupling with (the threshold condition) positive average 

gain at the bandgap edges corresponding to Re[Ake^/K} = * Tlj,/8, 

1 

The normalized temporal growth rate is Im{Ake ^/K] = - e./4e and 
' .a r 1 r 

the longitudinal wavenumber p = is real. (Bandgap edges will refer 

to the lossless case for this chapter.) 

Figure 5,4 shows similar mappings of the dispersion equation 

roots for gain coupling. In Figure 5,4a, b, c the merging root be- 

1 

havior occurs at the value Ak e ^/K = -ie./4e +1111. 1/8. Thus, for 

gain modiilation, absolute instability occurs at the center of the band- 
gap for average gain (e. < 0) and average loss (e^ > 0) whenever the 
threshold l'Hjl/2 is satisfied. The result for zero average 

gain is not surprising since the inverted bandgaps are similar to 
those of parametric instabilities and backward wave oscillators. 

However, note that gain -coupled media may also have instabilities 
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Fig. 5, 3 Locus of the roots of the dispersion relation in the Ap/K 

plane at the first Bragg order for index coupling with average gain 

a) e. /c = 0, 0, b) e. /c = 0, 03 and c) c, /e = -0, 03. Locus of 
^ r 1 r 1 r 

Ak e^/K for the above is given in d), .Similar loci in the Ap/K plane 
are produced for Re[Ake^/K] < 0. 








1 « /( = 0.0 
I r 


I t I 


€./c =0.05 



7= i 0. 1 


1 '“' 

^ 0 01 



O.Ol 0.02 


Re|A0/K| 


Re(Akv/<; Al 


Fig. 5.4 Locus of the roots of the dispersion relation in the 
Ap/K plane at the first Bragg order for gain coupling with 
average gain a) e.. /e =0, 0, b) e. /c ~ 0, 05 and c) c. / £ = “0, 05, 
Locus of Ak c ^/K for the above is given in d). 





when the average gain is negative as long as the gain periodically takes 
on positive values (i. e. the average loss i-s less than some fraction of 
the variation). 

As evidenced by the structure of the dispersion relations of 

chapters III and IV and the repeated Brillouin diagram of Figure 5. 5, 

index and gain coupling at the second Bragg resonance display similar 

characteristics, to the first-order index coupling diagrams (i. e. no 

inverted bandgaps appear). Thus, the stability of both index and gain 

coupling atthe second Bragg resonance should be similar to thatof the first 

order index coupling case since the stability criteria is dependent only upon 

the dispersion relation. The only difference should be in the sign of the 

bandgap shif^and hence the relative instability frequency which is positive 

or negative for index or gain coupling respectively. Indeed this is true. 

The mapping of the complex roots (not shown) of D(Ap, co) as w.-* 0 resembles 

Figure 5.3. The absolute instability occurs at the bandgap edges for index 

1 . 

or gain coupling with temporal growth rate given by Im[Ake^/K] at the 
bandgap. The threshold condition requires positive gain (e.<0) for abso- 
lute instabilities at the second Bragg resonance for both coupling types. 

The results of odd and even Bragg resonances show different 
stability characteristics as exemplified by the specific cases of N= 1,2 
discussed above. The mathematical conditions for absolute instability 
parameters can be found from the above considerations and the ECW 
dispersion relation. For any Bragg order N, oscillation takes place at 

Re{6j^/K] = ± Re{Xj^/K} (5.B1) 

with temporal growth rate 

Imi 6j^/K} = + Imfxj^/K} (5. B2) 

at or above the threshold 
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Im{Akj^^2/K} = 0 . (5. B3) 

The results for the first few Bragg resonances are summarized in 
Table 5. 1. It is noted that with increasing Bragg order, the 
threshold approaches zero for all couplings since r| is usually small. 
Also the instability frequencies tend to cluster above (index coupling) 
or below (gain coupling) the exact Bragg resonance for higher Bragg 
orders. This is due to the bandgap shift. 

The first two Bragg resonances have been shown to be arche- 
typical of all Bragg resonances. In particular we found that for index 
coupling, absolute instabilities occur for positive average gain at 
both bandgap edges with temporal growth rates given by co^ = o at the 
bandgap. For gain coupling, the behavior is identical to the index 
coupling case for N even, with the exception of the instability fre- 
quency which is shifted due to the bandgap shift. For gain coupling 
with odd N, absolute instabilities occur only at the center of the band- 
gap with temporal growth rates equal to the value of CO^ at the bandgap. 

The correct choice for the sign of Im[Ap/K3 depends upon the 
position of the roots of the dispersion relation as “* 0 or at the 
merging root point = cr. If the roots have crossed the real Ap /K 
axis the outgoing waves are amplifying and if they do not cross, the 
outgoing waves are decaying. Thus, for the case of absolute instabil- 
ities, there will be both amplifying and decaying waves at frequencies 
adjacent to the instability (see Figure 5. 3-4). However, this behavior 
is overshadowed by the absolute instability and is not important for 

infinite media. The proper sign of Im{6/K] is chosen by the fact that 


Z 


Coupling Order Threshold I Frequency I Growth Rate 




































the dominant time behavior is controlled by the highest root in the 
CO plane. This is equal to the value CO^ = o- that occurs when the 
roots of Ap/K merge. 

I 

In all cases where the gain is not large enough to support 
absolute instabilities, the medium is stable and the outgoing waves 
decay spatially due to coupling and losses. Thus, because there are 
no axis crossings of the roots of D(Ap, iO ) in the Ap/K plane, the sign 
of Im{Ap/K} > 0 is specified for outgoing waves. The sign of 
Im{6/K} is cho sen to correspond to the least lossy wave (i. e, 

Im{6/K} > 0) since this root produces the dominant field contribution. 
Thus, for stable media the periodicity will enhance the spatial decay 
and not affect the temporal decay for index and even-order gain 
coupling. For odd-order gain coupling the spatial and temporal 
decay are both diminished. 

C, Complex Coupling and Multiharmonic Periodicities 

In practical cases, both index and gain periodicity will occur 

together ('H = "nj. "nj)* This case, referred to as complex coupling, 

will usually cause absolute instabilities to occur at two frequencies 

due to finite- hj.* The explicit conditions are involved for arbitrary 

Bragg order but reduce to the following form at first order: absolute 

instabilities will occur in pairs at the bandgap edges (Re {Ak, e ^/K 3 = t "H /8) 

1 r X 

i 

with a growth rate Im.fAkj^ = -C./4e^+ the threshold is 

e./e^ = and thus cannot be reduced below the gain coupling value. 

At even Bragg resonances, where the threshold was pre- 
viously zero, a reduction in threshold can take place for the case 
of complex couplings. As an example consider the second Bragg 


resonance. Figure 5.6 repeats the Brillotiin diagram for complex 
coupling and = 0* For equal amounts of index and gain 

coupling (i, e. T] = T].) the curves give inverted bandgaps. Thus, 
one expects absolute instabilities at exact Bragg resonance. It 
can be shown that the instability threshold is reduced to 
e./e^ = 5 |t|^T|^|/ 6 (average loss). Two absolute instabilities occur 
for T) 4 and e./e =0 at the frequencies where Im{A(3 /k} is 

IT 1 1 IT M 

zero and where it is maximum. Thus, as the ratio approaches 

zero or infinity, the instability frequencies merge to the center of 

the bandgap apd as -* 1 the instability frequencies tend to 

merge in pairs toward exact Bragg resonance. The growth rates 

are found from the value of ImfAk^ s^^/K} at the bandgap center. 

The application of the stability criteria to multiharmonic 

periodicities is similar to the previous analysis (chapters III, IV) 

with the proper phase mismatch and coupling Xjq substituted into 

the dispersion relation. Although there is little or no effect at the 

first Bragg order for changes from sinusoidal periodicities to other 

typical periodicities, the effects at higher Bragg orders (N ^ 2) may 

be significant due to Fourier components f^^. It can be shown 

2 

(Appendix B) that the coupling is increased from values 0(“n ) for 

sinusoidal periodicities to values ©(T)) for sawtooth periodicities 

at the second Bragg resonance. 

At the third Bragg order, the coupling can be increased 
3 

from values 0(r| ) to values 0(ri) for square-, triangular - or 
sawtooth -wave periodicities. This implies that instability 



Fig. 5.6 Brilloiiin diagram at second Bragg order for ECW theory 
when both index and gain/loss coupling are present. 
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thresholds can be changed at higher Bragg orders by the use of non- 

II -2 - 5 

sinusoidal periodicities. Typically, |T]p~O(10 -10 ) in integrated 

optics applications. 


D. Comments on Stability 

In this chapter we have found the stability criteria for time- 
independent periodic media. This media may support either decay- 
ing (stable) or temporally growing (absolutely unstable) waves in 
basically dispersionless dielectrics. The ECW equations provided 
explicit values of threshold, instability frequency and temporal 
growth rates for absolute instabilities. Average positive gain was 
required for oscillation at the bandgap edges for index coupling and 
even-order gam coupling. For odd-order gain coupling, the average 
gain could be negative if the gain periodically took on positive values 
and the oscillation took place at the bandgap center. In some cases 
both types of coupling or multiharmonic periodicities could reduce 
thresholds. 

There has been some question as to the effect of boxindaries 
upon DFB oscillation. It is apparent from the preceding sections 
that instability takes place when certain modes of the active peri- 
odic media achieve threshold. No boundaries are needed since 
feedback action is produced by the periodicity. 

It appears that for frequencies other than the instability 
frequencies, waves may either decay or be convectively amplified 
if the structure is short enough such that the absolute instabilities 
do not occur. In this case, the active periodic structure might be 


X 


used as a filter -amplifier. However, each individual system must 

be considered along with the boundary conditions to determine sta- 
35 , 73 , 74 ' 
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The results of this chapter have been summarized. 
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CHAPTER VI 


APPLICATIONS OF ECW THEORY 

The results of previous chapters have depended only upon 
the dispersion diagram of xinbounded periodic media. In this chapter 
several cases of bounded media will be discussed. The use of the 
ECW equations allows behavior at all Bragg orders to be approxi- 
mated easily. The purpose of this chapter is not an exhaustive 
coverage of periodic structures but rather an indication of a wide 
range of problems which may be solved by the use of the previous 
theory. 

In section A, longitudinally bovinded passive media (i. e. DFB 
filters) will be covered along with examples of transients in a periodic 
slab. In section B the characteristics of higher-order DFB lasers 
will be given. Section C will cover the case of holographic grating 
diffraction and the last section will characterize the propagation of 
Gaussian beams in periodic media. 


A. DFB Filters 

1. Effect of Longitudinal Bovmdaries 

We consider waves in a periodic slab of length JL, The ECW 
equations account for Bragg coupling or scattering from the perio- 
dicities but not for the coupling or scattering due to the boundaries. 
The equations can be easily modified. The relative dielectric con- 
stant is now 


€(z) = e + ie. + rect(z/£) e Zj (r| +iri.) f cospKz 
r 1 r p r 1 p p 




where 
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rect(z/i) 



\z\ <1/2 

\i.\> 1/2 


For simplicity, consider singly periodic media with no average loss 


or gain (i. e. f ^ = 0 for p ^ 1, = 0). The additional Fourier com- 

ponents available for coupling or scattering arise from rect(z/2) and 
have not been accounted for previously. These components are pro- 
portional to the Fourier transform, 

00 

T] J rect(z/i) cos Kz e^^^ dz 


■00 


(6.A1) 


r\i /sin(0+K)£/2 , sin{p-K)i/2 1 
= 2 t“(^)t72~ (p-K)i/2 ■ ; 


Thus, the additional available perturbation per unit length is propor- 
tional to the right-hand side of (6.A1) divided by the length As in 
the case of multiharmonic periodicities, only the strongest effect 
(i, e. first order in Ti) of the boundary perturbation of (6. Al) upon Fj 
and Bj are considered. For higher»-order Bragg interactions, the 

coupling diagram of Figure 6.1 is helpful. Not only is there the usual 

N 

cros s- coupling through terms 0(ri ) and self-coupling through terms 
2 

0(T1 ), but there is cross- and self-coupling proportional to and T)0, 
respectively, where 


$ 


© 


. sin((3+K/2)jt sin(3-K/2)Ji 

(p+K/2)je (p-K/2)Ji 

, ( sin{p+K)f/2 sin(P-K)f/2\ 

I (p+K)jt/2 (p-K)f/2 J 

J 


(6.A2) 


2 

are the boundary effects. This is due to the terms [-k e^Ti(©Fj+ IBj)/2] 
and ( 0 + # F^)/2] which are added respectively to the Fj 

and Bj equations of the ECW equations (3.A6), Consequently, the 
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modified ECW equations become identical to the former ECW equations 
with the chcinges 6^^ -« and ”* Xj^j where 


^N _ ^N- _ N^ 0 ^ 

K ■ K ^ 






^N 


K 


+ f 

K 8 * 


(6. A3) 


In each case the first term of the above equation expresses the effect 
of the periodicity while the second term expresses the effect of the 


boxindary. For higher-order interactions, the latter effect predomi- 


,1-N 


nates in coupling when KX « (t)/2)'‘ while the opposite is true for 
,1-N 


Ki » (t)/2) This is expected since for thick slabs the extra 

Fourier components introduced by the trvincated periodic media are 
tightly clustered arovmd the components introduced by the periodicity 
of the infinite media and no new effects are observed. However, for 
thin slabs, the truncation produces Fourier components at many 
multiples of the Bragg wavenumber which are capable of directly 
coupling Fj to Bj. In this case the boundary coupling (« r|§) predomi- 
nates. Note that the above inequalities denote the region where the 
results of the infinite media can be directly applied to the longitudi- 
nally bounded case. Also note that the boundary effects can be reduced 
if the perturbation is gradually truncated at the slab ends. A similar 


effect has been noted in quantum mechanical scattering. 


79 


The 


boundary effect seldom affects the phase mismatch since it is negli- 

-1 


gible when KX » T| 


2. ECW Reflection and Transmission Coefficients 


In this subsection we consider the'general case of transversely 
and longitudinally bounded structures where each coupled wave repre- 
sents a different mode. Hence, the longitudinal wavenumbers may be 

3 1 

different for each mode. For simplicity the numerical examples 

will correspond to the special case of coupling between waves of 

identical modes in thick slabs although the analytic resvdts are given 

for tlie general case using the ECW approximations. 

Figure 6. 2 schematically shows the configuration. Consider 

ippZ 

a positive phase velocity wave F,(z) U (x)e ^ which couples to a 
negative phase velocity wave Bj(z)U^(x)e ^ . The subscripts p, q 
refer to the mode number of the transverse distribution U(x). The 
phase match condition for significant coupling at the Bragg order 
is 


Pp + Pq-“NK (6.A4) 

« 

where K = 2 it/A and A is the f\mdamental spatial period. The 
boundary conditions are found from Figure 6. 3 by use of the ECW 
equations which correspond to continuity of the electric field at 
z = ± J^/2. 


^-ipX/2 


+ 




ipjl/2 


n=-l l-2n/N 


+B (z=-X/2) e 
l-2n/N 


i(l-2n/N)i/2 


(6.A5) 



INCIDENT PULSE 


REFLECTED PULSE 


-i2/2 



PERIODIC STRUCTURE 

^ 


TRANSMITTED 

PULSE 


SL/1 


Fig. 6. 2 Configuration of pulses and periodic slab 
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N n^.l l-2n/N 

+ B(z=i/2)e-'(l-2-/N)V2 

•l-2n/N 


where 


N /2 = 



+ (J) S(z=i/2) 



for N even 

for N odd 
for N = 1 

for N even 

for N odd 


(6.A6) 


r 


V 


and and are the reflection and transmission coefficients. 
Similar equations result from matching the derivatives (or the mag- 
netic field) at z = ± Jl/2. The only additional equations needed are the 
ECW equations modified for transversely boxinded media. 


-Bj(z) - ir^BjCz) = iXqp 


Bj(z) 


Fj(z) 


(6.A7) 


where the coupling and phase mismatch are assumed to be known. 
Instead of solving the above equations exactly, a perturbation scheme 
is introduced which significantly simplifies the calculation and elimi- 
nates the bovindary coupling which has been accounted for previously. 

Introduce the ordering parameter X which will eventually be set to 
77 

unity. Consider a power series solution of the form 

F 


1 


= + Xf|^^ + X^ +. . . 


B 


1 




1 


1 


T 


+ . . . 


(6.A8) 
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€(Z) 


p oi(l-2/N)/9oZ 

n-2/N e 


F| 


Rn e 


B, e ‘V 


z = -£/2 


Fig. 6.3 ECW waves in periodic slab for coupling 
between identical modes (i. e. 3 =6 = 3). 

p q 
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for the waves and let the smallness be specified by \r). K all the 
waves are solved in terms of F^ and Bj (see equation (3. A6) ), then 
the electric field continuity equation becomes. 


•ipX/2 ^ ^i^in e'"P°‘^'^[l+0(\Ti)] 

^ ^ J|z=-il/2 


(6.A10) 




ip ^/2 ^ [f{°^+Xf|^4x^F^j^V. . .] le"^o^^^[l+0(XTi)] 

|z=£/2 (6. All) 

* *] le'^*^o^'^^[l+0(XTi)] 

'z=f/2 


[b|°4xb|^^+x^b^ 


where 


P = Po-^^^N 


p^ = NK/2 


,th 


near the N Bragg resonance. By equating terms 0(X) and phase 

19 

matching, the approximate bovmdary conditions become 


Fj(-jC/2) = 1 


Bj(-f/2) = 


Fi(V2) = 


Bj(je/2) = 0 


(6.A12) 


The reflection and transmission coefficients are solved using the 
boiandary conditions and the ECW equations. The details are carried 
out in Appendix C. We find 


R 




N 


N 


SLP_ 


Dj^cothDj^ 


.-i( 


-^r.+ b 


(6.A13) 


2 J 




Dj^cosh Dj^X- i^—Bl^)sinhDj^f 


(6.A14) 


where 


D 


N 


j- ,6 s-,- 


If the two coupled modes are identical, then 6 = 6^ = ^ and 

P q 

Xpq = Xqp = X* The subscript N has been dropped on 5 and x 
convenience. 

From the above equations for and Tj^ we observe the 
following ; 

1. maximum reflection 


= i tanh(x X ^ I 
TTjmax pq ^qP 


(6.A15) 


2. minimum transmission 
T 


N min 


= l/cosh(x X ^ > 
^pq^qp 


(6. A16) 


3. = 0 and Tj^ = 1 for 

V^q " ^^^pq (n=l,2,...) (6.A17) 

where the corresponding phase of is ( 2 n-l) 7 r/ 2 ; 

2 2 

4. IRj^I + iTj^l = 1 for pa ssive and lossless media. 

These results are demonstrated in Figure 6.4 for the case \ f = 2 
where the bovindary effects are neglected. The magnitudes of the 
reflection and transmission coefficients are shown as a function of 
6X(i. e. frequency deviation from bandgap center). Note that the 
half-width of the main reflection maximum is equal to the bandgap 
width. Hence, there is significant reflection outside of the bandgap 
due to the deviation of Re[Ap/K] as well as Im{Ap/K} away from their 
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tinperturbed values. Figure 6. 5 shows the eqviimagnitude and equi- 

phase curves of the reflection coefficient plotted as a fvinction of 

and '6^ where again the boundary effects are neglected. Note that 

Figure 6. 5a correspoiids to a horizontal line = 2) across Figure 

6. 4. If the coupled modes are different (i. e. p q), the curves in 

Figure 6.' 4 and 6. 5 are valid if 6 is replaced by (6 +6 )/2 and y hy 

i P ^ 

(y?^ yqp)^. Similar results for media with average gain have been 

discussed for the first-order index coupling case when the modes are 

identical. 

It is of interest to note the change in the reflection coefficient 
with Bragg order. For small coupling y^^J^ « 1 , 




(6.A18) 


Thus, the variation of with is 

RnI “ i{(Tl/2)^Kj^+Ti sin[(N±l)K£/2]} (6.A19) 

I max 

where the first term expresses the dependence upon the periodicity 

perturbation and the second term expresses the effect of the boundary 

\ 

perturbation. 

For other typical periodicities, the results of Appendix B will 
provide appropriate coupling or phase mismatch. In particular, note 
that for thick slabs, the coupling is constant at odd Bragg orders for 
square-wave periodicities. This implies that the maximum reflection 
coefficient is also constant at odd Bragg orders. These results are 
in agreement with well-known exact results that use matrix calcula- 
tions. Figure 6.6 exemplifie's th*ese concepts for three typical 


V 




K^/2 


K£ 

V .-/1 £ 


3K£/2 


i 


Fig. 6. 6 Sketches of reflection coefficient amplitude for ~ 0. 5 
for a) sinusoidal, b) square-wave and c) sawtooth- wave perio- 
dicities. In all cases the boundary effects are neglected. Note 
small reflection at N Bragg order when f^^ = 0. The bandgap 
shift is too small to observe on this scale. 
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periodicities with index coupling. Note that periodicities with large 
Fourier components will have large coupling and reflection co- 
efficients at the Bragg order. Thus, periodicities with odd 
Fourier components (e. g. square- or triangvilar- wave) will have 
large reflection coefficients only at odd Bragg orders. 

The use of dynamically generated (e.g. acoustic waves or 
electro-optic effect) mioltiharmonic periodicities opens up the possi- 
bility of controlling the feedback strength in DFB filters. This, in 
turn, varies the passband. From the results of chapter IV it is ap- 
parent that the passband could vary as sin(0/2) where 9 is the relative 
phase between Fourier components of the periodicity at the second 
Bragg resonance (see Figure 4,5). If the botindary effects are sig- 
nificant or predominant, then the previous results will be modified. 

That is, for thin slabs, the reflection will be significantly increased 

1 

at higher Bragg orders when {risin[(N±l)K + Jj^/2] is no longer 

negligible. 


3. Born Approximation Reflection Coefficient 

In order to summarize the properties of reflections from both 

thick and thin slabs, it is usefvil to find the Born approximation to the 

reflection coefficient. Consider the wave equation for some transverse 

component of the electric field E(z). 

2 

+ k^e 1e(z) = -XTik^e cosKzE(z)=-XS(z)E(z) (6.A20) 

Ldz'^ rj r 

where the ordering parameter X has again been used to show the 
smallness of r|. Assuming a solution 

E(z) = E^°^(z) + X E^^^z) + X^ E^^\z) + . . . 
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we find 


00 

= J S(z') g(z.z') E^”Mz')dz' 


(6.A21) 


-00 


where 


E<®>(z) = e'P' 


g(z,z') = 


2(3 


and where p = k for this computation. A straightforward calculation 


shows. 




,( 1 ) 


Thus E' ' represents a wave traveling to the left which is produced 
by the wave E^^^ traveling to the right. The approximation which uses 
the terms E^^^ and E^^^ is known as the Born approximation. The 
Born reflection coefficient is then the coefficient of e"^^^ in E^^^. or 


R 


(1) _ iripf r sin(3+K/2)f sin(p- 


4 C (P+K/2“)A ■ (p 


/n iri k e 2 i 

p_(l) _ _r — $ 


-K/2)X J 


(6.A23) 


which is valid for ( T] k f'|/4 « 1 at the first Bragg resonance 
or for |ti/4|« 1 at higher Bragg resonances. 

} 

At the first Bragg order, for 1 « 1 and hi « 1, 
the Born approximation and the modified ECW approximation 
become equal. Explicitly we find 


R 


( 1 ) 


i T] k e 2 Jt 


R, 


(6.A24) 


max 


' max 
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at the bandgap center for Ki ^ 10. For higher Bragg orders, the Born 

I 

approximation produces a series in powers of r| which accovints only 
'for the depth of the perturbation (i. e. the end ’effects) and not the 
length over which it acts. A series in r\Ki is needed to account for 
significant periodic effects. This latter series occurs through the 
modified ECW equations. Thus, the Born approximation is useful 
at the first Bragg order or at higher orders where the slab satisfies 
the relation K£ « For thick slabs, the ECW theory provides the 

proper phase mismatch, . coupling and reflection. 

4. Transients in Periodic Slabs 

Consider an incident pulse f(t) with frequency spectrum 

00 i . 

F(co) = J* f(t)e "Mt (6.A25) 

-00 

where the reflected and transmitted ptilses are, respectively, 

r(t) j”F(aj)RN( 0 ))e-^‘^ dw (6.A26) 

' -00 

t(t) = ^ r F(w)Tj^(a3)e''‘*^^d0J (6.A27) 

-00 

This case was shown in Figure 6.2. Numerical inversion has been 
used by a number of authors to obtain time response to radiation and 
scattering problems. The Cooley-Tukey fast Fourier transform 
(FFT)^^ has been used in the numerical examples with 2^^ = 2048 
samples to calcvilate r(t) and t(t). 

The cases of reflection and transmission of rectangular and 
Gaussian pulses of several center frequencies and widths have been 
carried out. Normalized pulse lengths T were chosen to be 0. 25 
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and 2. 0 time vinits where each time unit corresponds to the transit 
time of the pulse across the slab length 1 . As in previous sections, 
the numerical results are valid at the first Bragg order and all higher 
Bragg orders where the boundary effects can be ignored. 

Figure 6. 7 displays the envelopes of the reflected and trans- 
mitted pulses which result from an incident pulse of length T = 0.25. 
Several values of coupling are used and the carrier frequency is at 
the bandgap center (6^^ = 0). In all the illustrations, the time value 
t = 0 corresponds to the instant when the center of the pulse is at the 
first boundary. For weak coupling, the reflected pulse is spread over 
2.25 time vinits (Figure 6.7 a,c) because the energy is reflected from 
successive striations and the echo from the last boundary has a rovind 
trip time of 2 units. Since the successive reflections are relatively 
weak, the reflected pulse is quasi- rectangular while the transmitted 
pvilse is similar to the time delayed incident pulse. As the coupling 
is increased, multiple interference leads to transients over a longer 
time period (Figure 6. 7 i, j). The reflected pulse reaches a maximum 
^ at 0. 25 vmits when the entire incident pulse has just entered the slab. 
The subsequent fall-off in the reflected pulse amplitude is due to the 
fact that a large portion of the signal has already been reflected due 
to the first few striations. For large coupling, the transmitted pulse 
consists mainly of two peaks, 0.25 units apart, which correspond to 
differentiation (i. e. high-pass filter) of the incident pulse. 

The same coupling strength sequence is given in Figure 6. 8 
for a rectangular pulse of length t = 2. For weak coupling (Figure 
6.8a), the reflected pulse is similar to the autocorrelation of the 
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REFLECTED OUTPUT TRANSMITTED OUTPUT 



NORMALIZED TIME 


Fig. 6. 7 Reflected and transmitted pulses for different values of 

coupling coefficient. Incident rectangular pulse has length 0.25 

time units and a carrier frequency at R--i . Note the difference 

' max ■ 

in vertical scales. 
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incident p\ilse. This indicates that the reflection coefficient is similar 
to the incident piilse spectrum. For strong coupling, the transmitted 
pulse is again similar to the absolute value of the derivative of the 
incident pulse with characteristic interfering echos. Note that the 
reflected pulse may have an amplitude larger than unity (Figure 6. 8g,i) 
due to constructive addition of successive reflections. 

Figure 6.9 displays the reflected and transmitted pvilses for 
a Gaussian incident pulse of width T = 2 (width is taken at the 1/e 
values). Since a Gaussian pulse contains smaller higher frequency 
spectral components, the reflected and transmitted pulses are grossly 
similar to the incident pulse. For narrowed Gaussian pulses (not 
shown here), the results are somewhat similar to those obtained for 
a rectangular pulse of the same width. 

Several examples of pvilses with carrier frequencies at the 

first several zeros of the reflection coefficient have been given else- 
29 

where. In general, the transmitted pulse is similar to the incident 
pulse. However, the reflected pulse is broken up into a long pulse 
train which can be considerably longer than the incident pulse width. 

5. Discussion of DFB Filters 

In the preceding subsections, the ECW theory was applied to 
longitudinally bounded media. The ECW equations were modified to 
directly account for the coupling between Fj(z) and Bj(z) due to the 
boundary. This changed the coupling so that it consisted of two parts 
with the following physical meaning. One part of the coupling, due 
only to the perturbation is proportional to (“n and hence increases ' 









I lh^r~nrn’iilivi^i*il II 
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with slab thickness since it is due to repeated reflections in the di- 
electric. The second part of the coupling is due to the above mentioned 

X 

boundary effects and is proportional to [t) sinC(N±l)K+Ake^ 3 This 

accovints only for the interference effects due to the periodicity trun- 
cation and may become zero for certain lengths at particular frequen- 

V 

cies. The application of the boundary conditions was then introduced 

in such a way that end effects were not accounted for again. 

The validity of this approach was tested by using the exact 
89 90 

results of Bedrosian. * The comparison was made on periodic 

media of length Ki> = 26ir with Tj = 0. 05. The ECW and exact invariant 

imbedding technique agreed to within 1^ in the two cases where the 

-2 

boxindary effect was insignificant (R, i = 2. 55x10 ) and where it 

^ Imax 

7 90 

had a large effect (Ro 1 “ 5.88x10"^) at the second Bragg order. 

'max 

Comparisons to the Born approximation showed that at higher 

N 

Bragg orders, the Born approximation ignored the T) KX’ term of the 
coupling. However, at the first Bragg order the Born and ECW ap- 
proximations are in very close agreement. For cases where Ti -* 1, 
the full Floquet theory must be used. This is straightforward in 
principle but laborious to carry out. The disadvantage is that lengthy 
numerical calculations have to be carried out and the intuitive appeal 
of the ECW theory is lost. 

We have shown that periodic slabs may be useful as edge 
differentiations or multiple pulse generators. The transient re- 
svilts may be useful in the fields of microwave and optical filters 
and for sounders of subsurface layers. 
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B. Higher-Order DFB Lasers 


54 


Ever since Kogelnik and Shank's demonstration and descrip- 
23 

tion of DFB lasers in the early 1970' s there has been great interest 

in this field. Since then much work has been done at the first Bragg 

order. in 1972 Bjorkholm and Shank 

demonstrated the first higher-order DFB laser with output at the 

second and third Bragg resonances. More recently there has been 

71 

further experimental work at the second Bragg resonance and with 

70 

multiple frequency DFB lasers. 

In this section the approximate mode spectrum cind threshold 

gain are given as a function of Bragg order. The results, which use 

the ECW equations, are analogous to the coupled wave results of 

23 

Kogelnik and Shank. Hence, the derivation will be only briefly 
sketched here. The details are given in Appendix D. 


1. ECW Result for Threshold Gains and Mode Spectra 
Oscillation will take place when output occurs for zero input. 
This can be .alternately stated as the condition where the reflection 
or transmission coefficient becomes infinite. We use the ECW re- 
flection coefficient discussed in previous sections and set the denom- 
inator equal to zero or 

Dj^ coth Dj^A = iA6^ (6. Bl) 


This can be put into the form 
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r°N"‘ S 


-1 


(6.B2) 


The complex solution to the above transcendental equation produces 
the threshold gain and the mode spectrum when solved for amplitude 
and phase. 


2. High -Gain Approximation 

The high-gain approximation assumes that-the periodicity 
has little effect on the propagation, or, for 


- 

K 




K 


32(N^-1) 


r ^ ^ 8(N -1)-* 


K - K “ K 


(6.B4) 


the condition is 


V « 
Xn 


g 


N 


(6.B5) 


This is approximated for singly periodic media by the inequality 


lr|/2l^ «I e./e^l (6.B6) 

From Appendix D, the threshold gain approximation is given by 

2 2 ’i' 

and the longitudinal mode spectrum is given by 

tan"^ ~ + phase(Xj^) m=0, ±1, ±2, . . . 

(6. B8) 
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where the first term is usually small and can be neglected. These 
expressions neglect boixndary effects which can be accoxinted for by 
a Straightforward modification as in the previous section. Equation 
(6.B8) is first solved for the normalized frequency deviation from 

the bandgap center and then (6.B7) is solved for g^^, the normalized 
average gain. It is apparent that for singly periodic media, the 
average gain will increase drastically with Bragg order. However, 
boundary effects or multiharmonic periodicities may reduce the 
threshold gain by increasing the coupling. 

The .mode spectrum is asymmetric and is sketched in Figure 
6. 10 for higher Bragg orders. For index coupling (Figure 6. 10a), the 
gain symmetrically pushes the modes outward from the usual two 
mirror cavity case (shown by dashed lines). However, the bandgap 
shift produces an asymmetry and shifts the entire spectrum toward 
higher wavenumbers. As in the first order case, no oscillation 

/ 

takes place at exact Bragg resonance. Figure 6. 10b, c shows the 
analogous information for gain coupling at odd and even Bragg orders, 
respectively. For odd Bragg orders, oscillation is possible near 
Bragg resonance, but the bandgap shift prevents oscillation exactly 
at A 1 = 0 , As expected from the Brillouin diagrams, even-order 

gain/loss coupling is similar to even-order index coupling except for 
the sign of the bandgap shift. 

3. Low-Gain Approximation 

Low gain implies the inequality for singly periodic media 

« |,Ti/2 1^ 


(6.B9) 


THRESHOLD GAIN 




-37t/2 


-7T/2 


7T/2 


37T/2 


AkN/c^i 


Fig. 6. 10 Sketch of mode spectrum in higher-order (N 2) DFB 
lasers for; a) index coupling (N = 2, 3, 4, . . , ); b) gain coupling 
(N = 3, 5, 7, ... ); and c) gain coupling (N = 2/4, 6, , . . ) in the high- 
gain approximation. 
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Again the details of the approximation are given in Appendix D. 
Boundary effects are neglected since lov? gain implies a large 
coupling -length product. Hence the boundary effects are small. 

Consider first .the case of index coupling. We find that 
the lowest -order longitudinal mode frequency is given by 



which becomes 




a. i 




(6. BIO) 


(6.B11) 


for long structures. Therefore, the oscillation takes place near the 
bandgap edges. This agrees with the stability analysis of the pre- 
vious chapter as expected. Identical reasoning holds for the case 
of gain/loss coupling at even Bragg orders. 

The threshold gain condition for index coupling and for even- 
order gain/loss coupling is given by 


g-IVT ^ 

&N 




(Xn 


(6.B12) 


This is approximated by 


e. 

1 


e 


1 

(K£)^ 


(6.B13) 


r 

for singly periodic media. Hence the threshold gain varies inversely 
with (ri/2) ^ for singly periodic media and inversely with the length 

cubed. Thus the threshold gain condition predicts that the average 
gain approaches zero as the length becomes infinite. Again this is 
in agreement with the stability criteria results. 
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For gain/loss coupling at odd Bragg orders the lowest- 


order longitudinal oscillation takes place' at the bandgap center 


(6.B14) 


with the threshold condition on the perturbation 




(6.B15) 


for zero average gain. Thus for DFB oscillation to start with = 0, 


the necessary critical length varies as 


I 

C 


(hl/2)^ 


(6.B16) 


with Bragg order. The mode spectrum (6,B14) agrees with the 
stability prediction and (6,B15) implies that negative average gain 
can produce oscillation oi l> JL^. This is verified by the approximate 
threshold gain condition on the average gain which is 


and becomes 


(6,B17) 


I 


(6,B18) 


as the length increases. For singly periodic media the threshold 
gain varies as (r|/2)^. ' This also agrees with the stability criteria 
predictions since this condition is identical to equation (5, B3). 
Oscillation again takes place at the bandgap center given by 


0 


(6.B19) 
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Note that only the lowest-orcler (m = 0) longitudinal mode 
characteristics are given in the low -gain approximation. 


4, Discussion of Higher-Order DFB Lasers 

In all cases, the previous results are extensions of the well- 

23 

known work of Kogelnik and Shank to higher Bragg orders. The 
ECW theory gives the correct and 6^^ to use in the theories. 
Numerical results have been given at the first Bragg order for both’ 
transversely unbounded^ and bounded^"^' media and so can 
be used at higher orders for guided modes. 

The ECW theory shows that the mode spectrum is asymmetri- 
cally shifted from the exact Bragg resonance for high gain and has 
characteristic differences which depend upon coupling type and Bragg 
order. Note that in the low-gain case, the threshold gain and fre- 
quency of the lowest -order longitudinal mode is also predicted by the 
stability criteria without regard for the boundary conditions. In this 
case, enough coupling or feedback is available to produce oscillation 
without boundary coupling. In general, threshold gains increase 
drastically with Bragg order unless multiharmonic periodicities 
are used or unless boundary perturbations play a significant role. 

The low-gain threshold varies as (11/2) or as (r)/2^for the case 
of index and even-order gain coupling or odd-order gain coupling, 
respectively^ for singly periodic media. 


C, Higher-Order Hologram Diffraction 


There has been interest in transversely varying periodic 
media (Figure 6.11). This configuration represents electromagnetic 
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Fig. 6. 11 Configuration of TE wave obliquely incident upon a 
holographic grating. The original wave and the Bragg reflected 
wave emerge as and respectively. 
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wave diffraction by acoustic waves or by holographic gratings. The 

22 

former problem was solved by Chu and Tamir who used both the 
Floquet and the coupled waves approach. The latter problem was 
dealt with exhaustively by Kogelnik^^ who developed a coupled waves 
approach that was valid near the first Bragg resonance. More 
recently variations of these problems have been investigated theo- 
retically and experimentally by several authors. 

In this section the ECW equations are discussed for trans- 
versely periodic media. The derivation is similar to that of chapter 
III and the details are given in Appendix E. The approximate boundary 
conditions are applied and the results are given for the strengths of 
the xmdiffracted and Bragg diffracted waves in holographic gratings. 

As in previous chapters, the ECW results will hold for all Bragg 
orders. 


1. ECW Equations for Transversely Periodic Media 


Consider solutions to the wave equation for the transverse 


electric (TE) field 


2 2 

r ^(x)1e(x, z) = 0 

^ 9z 9x 


(6. Cl) 


where 


e(x) = e + ie. + e Z/{ti +iri.) f cos pKz 
' ' r 1 r p T ip p ^ 


is the periodic relative dielectric constant. Following the example 
of Kogelnik^^ we consider a wave of form which is suc- 


cessively scattered from the periodic dielectric to a wave of form 

-NK) x+k^ first wave is designated as Fj^(z) and the Bragg 


I 
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scattered wave as Bj(z). Using the same ECW assumptions as in 
chapter III we let 


n*/e, 

E(x,z) = S F(z)j ^3 ^ 
n=-l' ( 


^i[(l-2n)/N)poX+Apx+kzQzJ 


^ ®<^>l-2n/N (6. C2) 

+ (0> S(^) 

be the assumed TE field. Substituting the above expression into the 
wave equation (6. Cl) for singly periodic lossless media, we find the 
following ECW equations (see Appendix E). 


^1 ~ ^ % ®1^^^ 
Bj(z) +i 6j^Bj(z) = i Fj(z) 


where 


N 


k^e 1 -C^T ( o r~ o 

' "N^2sm0^- 

N 




2 k 


(6.C3) 


(6.C4) 




(-1)^^^ k"e 
2^+1 k, (sin^e 


*[4n(n-N)/N^}^ 


(6. C5) 


k 

P 

k 

2 

P. 


k + Ak- , 
o N 


P + ^Ptvt 

*^o 


k^ + Ak 
o 


NK/2 


2 2 2 
k e = kz +p 
o o 


sin 0 = p /k 

o *^o o 
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iAk ^ 

Assuming a solution of the form e ^ gives the dispersion relation 




(6.C6) 


The similarities and differences between the longitudinally periodic 
(LP) and the transversely periodic (TP) case are given below. 

1. The sign differences in the TP case cause the dispersion 
relation to be real for real periodicities and passive 
lossless dielectrics (i. e. = 0 = e.). Hence no bandgaps 
appear in the longitudinal wavenumber as in the LP case. 

2. Maximum phase matching between Fj(z) and Bj^(z) occurs 
slightly away from the exact Bragg condition for higher - 
order Bragg interactions in both the TP and LP cases. 

The first Bragg order results are found from (6.C4-5) as 



Ak j Vc 
K 


A6, sin 9 
^1 o 

K 


1 

COS0 

o 




-- “n 

8 cos 0 sin 0 


(6.C7) 


(6.C8) 


o o 

Small changes in frequency or angle are accounted for by the first 
and second terms of the phase mismatch. The third term is negligible 
away from normal incidence. Results for other Bragg orders and 
for complex periodicities can be found easily from previous results. 

A similar derivation for TM waves can also be’ made. 


V 


2. ECW Reflection and Transmission Coefficients 
The slab configuration of Fig. 6. 11 shows the undiffracted 
wave proportional to Tj^ and the Bragg diffracted wave proportional 
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to Rjyj» The incident wave is of unit amplitude. The appropriate 

' 15 

boundary conditions for the TP slab are approximately 


FM/?) = 1 
B(-X/2) = 0 


F(X/2) = Tj 
B(2/2) = R, 


(6.C9) 


These equations indicate that the Bragg diffracted wave increases 
in value from zero to while the undiffracted wave decreases from 
unity to Tj^. Note that any reflected wave in the region z < - X/2 is 
ignored. 

The ECW equations (6, C3) are solved with the boundary con- 
ditions imposed by (6, C9). The method of solution is analogous to 
that of the LP case (see Appendix C). The results are 

i X sin Ak„ H 
®-N = — Tk 

cos Akz-T sin Ak Z 

rr _ N // 


for coupling between identical modes (i. e. 6 = 6 , X = X = X)» 

Note that for exact phase matching, the reflection and transmission 
coefficients vary as sines and cosines. Hence, the energy is alter- 
nately shifted between the undiffracted wave Fj^(z) and the Bragg 
diffracted wave Bj^(z), along the z coordinate. Away from^'exact 
phase match the energy transfer is incomplete. 

3. Discussion of Holographic Gratings 

We note that the development of the ECW equations for the 

TP case is identical to that of the LP case. The ECW equations are 

' 19 22 69 

of the same form as that of several other authors, ’ ’ The 
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advantage of the present derivation is the simplicity of the explicit 
expressions for phase mismatch and coupling. In particular the 
deviation from the exact Bragg condition that provides maximum 
coupling or highest diffraction efficiency is easily found. 

In the preceding derivation we found coupling between the 
waves Fj(z) and This means that the diffracted and undiffracted 

waves, the amplitudes of which are given by and T^^, are sym- 
metric with respect to the z axis. A similar theory could be devel- 
oped for coupling between any other set of two waves. For example, 
one could find a set of coupled equations for coupling between Fj(z) 
and B(z)j^ 2 q/j^ where q = 0, 1, 2, , , . represents different spectral 
orders. In this case the vindiffracted wave propagates at the angle 

9 ^arcsin (P /k ) whereas the diffracted wave propagates at the 
^o 

angle 9 arcsin[p^(l+2q/N)/kz^] with respect to the z axis. The 
diffraction efficiency DE is usually defined by the relation DE = ITj^I . 

The theory is approximate since only a few space harmonics 
were used, only the approximate boundary conditions were applied 
and because the (assumed small) reflected wave for z < -£/2 was 
ignored. 


D, Gaussian Beams in Periodic Media 

With the advent of the laser, it is useful to consider Gaussian 

85 91-94 

beam scattering from periodic structures. Several authors ’ 
have considered Gaussian beam scattering at dielectric interfaces 
or in periodic media. The method of beam reflection at interfaces 

91 

was developed by Br ckhovskikh. We. will use this formulation and 


V 
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a similar one used by Tamir and Bertoni, 

In this section, the formtilation for Gaussian beam propaga- 
tion in LP and TP slabs will be outlined. The previously derived 
ECW reflection and transmission coefficients will be used. Numerical 
examples will be given in a future report. However, the general 
characteristics may be determined from the previous transient anal- 
ysis for certain cases. 

Consider the aperture Gaussian beam E^^(x, -h) of Figure 6.12 
that is formed at the plane z = -h. 


E (x>, 0) = 
ap' ’ ' 


-(x«/w/ 


w 


(6.D1) 


The propagating wave is given approximately by 


E^p(x, z) 


^-(xcos 0/w) ^ik a/F” [x sin9+(z+h)cos0] 


z=-h 


w 


(6. D2) 


' z = -h 


The incident wave i® given by the Fourier transform of the 

spectral amplitude $ (k^) 


E. (X. z) = f 5(k 
inc' , ' X 

2tr -00 


(6. D3) 


- h < z < 0 


where 


2 2 - 

k = (ke-k)^. 

z X 


The spectral amplitude is given in terms of the aperture field by 
the relation 


00 

$(k^) = ; E (X, -h)e“"‘^ dx 


(6. D4) 



Fig. 6. IZ Gaussian beam incident upon L.P or TP 
slab near Bragg resonance. 
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$(k ) = 

X 


^-[w(k^-k - n/c” sin0Q)/2 cos0q]‘ 
COS0 


(6.D5) 


Two assumptions are made. First, assume that the beam 
is well -collimated, or 


k//e~ w » 1 


(6. D6) 


where w represents the beam width. This limits significant spectral 
amplitude to values about k^ = k sin0^. Hence, k^ is usually real 
and the contribution to surface or lateral waves can be ignored. 

This is also necessary for the ECW approximations. Second, the 
incident and scattered waves are near Bragg resonance, or 


i i 

p /k £2 = N K/(2k e2) 
'^o o ' o 


f COS0Q 
I sin0Q 


LP media 
TP media 


(6.D7) 


The transmitted and reflected fields E„(x, z) and E„ (x, z) are 

J. Ix 

found from the integral of the product of the spectral amplitude 
(at z = 0) and the appropriate coefficient; 


E^(x, z) = ^ J°°T(k^)$(k^)e'^'^^'^^z^^'^^^^dk^ (6.D8) 

-00 

Ej^(x, z) = ^ dk^ (6.D9) 

-00 

where the ( signs hold for ( ^ p ) media. 

The expressions for R(k ) and T(k ) can be found from equa- 
tions (6.A13, 14) and (6. CIO, 11) for LP and TP media, respectively. 
The values of phase mismatch and coupling have been derived for 
arbitrary angle of incidence in TP media and at normal incidence for 
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LiP media. The extension of the LP results to arbitrary angles of 
incidence is a straightforward derivation and will not be given here. 


The explicit results near the first Bragg resonance are, 

2 


5 - Ak - tj £ Ak — I 

X “ [-ir- - TT 


1_ 
cost 


/^k' 
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for L.P and 
TP media 


(6. DIO) 


K 8cos0 f cosQ \ 

° Uin0°; 
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sin0 = 
o 


Ak 


c 

xo o ^ 

6 / k e2 
^o o 


for 
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for LP me^ia 
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o 

O' ? i 

[k'^e - p 
o ^o 


for LP media 
for TP media 


(6. Dll) 


K/2 


Note that for LP media the expressions for the transmitted 
and reflected fields are similar- to the transient analysis if the time 
coordinate is replaced by the space coordinate. In particular. Figures 
6. 7-9 show the spatial dispersion of pulses in space for rectangular 
and Gaussian input beams at exact Bragg frequencies and non-normal 
incidence if the normalized time is replaced by normalized distance. 

In this section the derivation for spatially bounded beams was 
given a form that is easy to compute under the ECW assumptions. 

It was shown that temporal and spatial dispersion are similar and 
that the results of the transient analysis could be extended to include 
beam propagation. 


- 153 - 


CHAPTER vn 

» 

CONCLUSIONS ' 

This report establishes an approxinaate method of calculating 
the properties of the Brillouin diagram at all Bragg orders for waves 

V 

in periodic media. The method is introduced in the second chapter 
where the connection between the Floquet and the coupled waves 
theory is shown and demonstrated with several numerical examples 
at the first Bragg resonance. In the third chapter,, the idea of cross- 
and self- coupling helps to extend the coupled waves theory to all 
Bragg resonances by the use of coupling diagrams. These results 
explicitly show the dependence of the bandgap width, bandgap shift 
and coupling upon Bragg order. The results closely match the exact 
Floquet dispersion relations. Furthermore, the ECW theory acco\ints 
for multiharmonic periodicities and demonstrates the idea of dis- 
appearing bandgaps under certain conditions. 

The fourth and fifth chapters deal with active or lossy dielec- 
trics. Inverted bandgaps occur only for certain types of coupling and 
at certain Bragg orders. The stability of active periodic media has 
characteristics that depended on the nature (i. e. inverting or non- 
inverting) of these bandgaps. Absolute stabilities are found to occur 
under certain conditions and only at certain frequencies in active 
periodic media. These threshold conditions and mode spectra agree 
with results fovind from an entirely different analysis, in the appro- 
priate limit, at the first Bragg order. The advantage of the stability 
analysis is that only the dispersion relation is needed to fully describe 
the stability characteristics. 
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The last chapter includes a few applications of the preceding 
theory to finite length media. The topics pf reflection, transmission, 
■transients, DFB lasers, holographic gratings -and beams in periodic 
media are briefly discussed and illustrated. 

The ECW theory demonstrates some of the power and versa- 
tility of the coupled waves formalism. It is anticipated that space- 
time periodic media and variable-frequency or almost-periodic media 

may be treated in a similar manner. The case of transversely 

' 32 

bounded media can be treated using a previously developed approach 
for the guided modes. The resiilt is that each power of rj will be mul- 
tiplied by an overlap integral. Hence the coupling will, be somewhat 
decreased. 

Exact theories such as the Floquet theory, matrix theory or 
the method of invariant imbedding will be useful for specific cases 
where exact results are needed. However, the intuitively appealing 
ECW theory gives results explicitly without lengthy computations 
and is surprisingly accurate in the cases treated here. 
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APPENDIX A 


ASYMPTOTIC FORM OF FIELDS FOR 
ABSOLUTE INSTABILITIES IN PERIODIC MEDIA 


Consider the following eqxiation for F(z,t) which is proportional 
to the electric field 


F^z t) - r“ r" d(A(3)dg; 

^ ' " J J D{A(3. ui) .2 

-00 - 00 '1-1/ ^2it) 


(A. 1) 


where S(Ap,co) is the normalized source and D(Ap,o>) is the dispersion 
relation. We are interested in the asymptotic value of F(z, t) for 
large times when an absolute instability occurs. This happens when 
the iToots of D(Ap,co) merge in the Ap plane from opposite half-planes 
separated by the line Im[Ap3 = 0. We assume that the source is 

analytic in Ap and that it is turned on at the time t = 0. 

73 74 

We follow the work of previous authors ’ and expand the 
dispersion relation about the instability at CO = CO' and Ap = Ap ' in a 
Taylor series. Then the integration is first carried out in the Ap 
plane and finally in the (0 plane. 

The Taylor expansion of D(Ap, co)/S{/ip, CO) is, 


where 


D(^P.CO) 

S(Ap,co) 


(^P - Ap')^ + D^ (CO-CO') 

i 9^(D/S) I 

3{Ap)^ |Ap=Ap' 


2 ^ 3(D/S) 

2 9co co=co' 


(A. 2) 


We rearrange the dispersion relation to find 
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F(Z, 03 ) = 


1 1 


4^iD^D^(03- 03')^ C 


i6pz 

e 

Ap-Ap'-iD^l W- aj')^/Dj 

(A. 3) 

iApz 


Ap-Ap' + iD2(0)- C0')2/Dj 


d(Ap) 


where 

F(z.t) = 




F(z, w) 


-itOt 

e 


2tt 


and the contour C is shown in Figure A. la. The first pole of (A. 3) is 
located in the upper half plane and contributes for z > 0 while the 
second pole is located in the lower half plane and contributes for z < 0. 
For z > 0, equation (A. 3) becomes 


1 

i[Ap'+i D2 (o}- w')2/D,]z 

F(z, w) = 1 (A. 4) 

2i Dj to- to ')2 


To carry out the integration in to, we depress the contour arovmd the 

singularities down to the real to .axis. This contour C" is shown in 

Figure A. lb where two poles at to = ± to represent the harmonic 

s 

source contributions. This second integral is 


F(z,t) = 


iAp'z 




21 D 1 D 2 


-[D2((0- tO')^/Di] z g-itOt 
: 

(t0-t0')2 


dcO 

2tt 


(A. 5) 


1, 

The branch cut of (to-to')^ provides the major contribution. Hence, we 
approximate the integral (A. 5) as t -♦ 00 by 
iAp'z 


F(z,t) = 


2iDi 


^ -[D2(<0-CO')2/Di] o +[D2(t0-t0')2/D]z 

/ ^ ^ J - 


(t0-t0')2 "o (to- t 0')2 

-itOt dtp 
2ir 


X e 


(A. 6) 



Fig. A. 1 a) Contour C in Ap-plane showing poles when 
Ap' = 0. b) Contour C" in a>-plane showing 2 poles due 
to source at xto and branch^^cut at instability frequency cO 
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i^P'z -iobt ^ cosh[D,(i0J.-ia)2/D. ]z idto. 

6*^6^ r 2 1 1 1 


F(z, t) = i D *f 

12 o 


(i 60^-ia)2 


2ir 


{A. 7) 


where £*)' = CO + io 
o 


Let (co. - a) = -q 


F(^,t) = J <=°°hCPl (-i5>^^°n^dq ,a. s) 

o (“i)2 q 2 


Ztt Dj 


Letting the upper limit tend to infinity we get 

i(A(3'z- cOot) ot 

F(z,t) - - — 7 V-j- as t-oo (A. 9) 

2-rr2 Dj D^C-i)^ t2 

Thus, the electric field is time-growing with the wavenumber and 
frequency of the absolute instability. For two instabilities at 
different frequencies, there will be two branch cuts m Figure A. lb 
and the F(z, t) will have two similar contributions, each of the form 
of (A. 9). 
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APPENDIX B 

ECW PARAMETERS FOR SQUARE- WAVE, 

TRIANGULAR- WAVE AND SAWTOOTH PERIODICITIES 

For applications to instabilities and DFB filters and oscillators 
it is convenient to have the phase mismatch and the coupling at 
the first three Bragg orders for square-wave, triangular -wave and 
sawtooth periodicities. 

th 

It has been shown that the major contribution to N order 
coupling is from the Fourier component f^^ of the periodicity. Hence 
for the three periodicities under consideration, the results at first 
and third orders are a straightforward application of the results of 
chapter III since each wave contains odd harmonics. In addition the 
sawtooth wave contains both even and odd harmonics so previous 
results can be used at all Bragg orders for this periodicity. The 
computation that needs to be carried out is for parameters at even- 
order resonances with odd-order Fourier components. We will carry 
out these computations for N = 2 in the lossless case (e = e^). 


The three periodicities have the following Fourier decompo- 
sition; 





p odd 


Square- wave 

f = ^ 
P 

1 

V. 0 

p even 

(B. 1) 

Triangular- wave 


^(-l)(p-l)/2^p2 

p odd 

(B.2) 



L 0 

p even 


Sawtooth 

f = 
p 

(-l)^P'^^) /p 


(B.3) 

with the normalization f. = 

1. An appropriate 

sine or 

cosine Fourier 


series is used for periodicities that are odd or even iri z, respec - 
tively. 

At the first Bragg resonance, we find that Xi/K = rj/8 and 
1 ' , ^ 

5j/K = tik. £^/K by the use of (3. A24, 3.A25) for all three periodicities 

This is identical to the sinusoidal case. 

At the second Bragg resonance, we find X^/K = -T) /8 and 
1 2 

6 /K = Ake^/K-ri /12 for sinusoidal periodicities. The other perio- 
dicities present need to be calculated. We consider the seventeen 
waves Bg, B^, . . . , S, . . . , F^, Fg for periodicities with fp = 0 for 
p even and N = 2. Thus, we need only the eleven waves Bg, Bg, . . . , 

S, . . . Fg, Fg since no significant coupling can take place through 
B^, F^ with q odd. Using the ECW approximations, we need to 
solve the following equations. 

-63 F, 


8 


-35 F, 


-15 Fv 


- 3 F- 


2Ak 


2 i Fj 

T 

k e 2 


2Ak 


- 3 B- 


2i Bj 

k 


= -^f F 
2 7^1 

+f3(F^+B^) +fj(F2+S)] 
= -£rf,(B^+Fj)] 

+ f3(B^+F2> +fj(B2+S)] 


(continued on next page) 
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15 B. 


-35 B, 




•63 B 


8 


H f B 
2 7 ^1 


(B. 4) 


2 

where self- and cross -coupling 0(r! ) has been included. The equations 
(B, 4) are solved to give the usual coupled equations for Fj and B^ with 


^ - ILf- f^ 

K 8 L 1 


—2 
°2 Ak e2 

K “ K “ 12 


2 ^^ 1^3 2 

+ — ^ + f. f. + 


ff 2 . i 

L 1 5 


.2 

‘3 


-i-f f +' 1 

35 5 7 •••J- 

(B. 5) 

2 


‘5 "210 • • : J 

(B.6) 


In (3. B6), if the terms drop off as f^ “ 1/p or faster, only the first two 

terms need to be kept in the series for accuracies of order 1^ in phase 

mismatch. In (3. B5) the series can be summed explicitly by using 

98 

Schlafli's Polynomials^ of order n for the square- and triangular- 
wave. Two terms will again give accuracies of order 1^. For square- 
waves we find 


K “ ' ^ ^2^ ’ K 


Ake 


^(.911) 


(B. 7) 


while for triangular waves 

X 


K 


where = 


2 6, Ake 2 


tt^-8 

“T5" 


32-3u 

■“ 64 " 


0. 1168 ... 

“ 0.3736 .... 


- ^(.990) 


(B. 8) 
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Thus, we have accounted for the major self- and cross- 
2 

coupling terms ©(t] ) which can be approximated by only considering 
the seven waves Bp S Fj, F^, F^ fdr periodicities with 

odd Fourier components f^. 

The case of the sawtooth- wave must be treated separately 
since it contains both even and odd Fourier components and must be 
expanded in a sine series. To consider the correction to the self- 
coupling from rif j and to the cross-coupling from we consider 
only the first three Fourier components fp f^, f^ and the waves 
B^, Bp S, Fj^, F^, F^ for simplicity. The equations for N = 2 

in the ECW approximation are then. 


-15 F4 


H f F 

2i 3 ^ 1 

1 

00 


= -,:?^f, F, 

2i 2 1 

-3 ^2 


= - f F 

2i 1 1 

2 ^1^ F,. + 
k ' 1 

21 f; 

k e ^ 



k 

-8 B3 

-15 B, 


^1 ^ 


2i F, 


= -i^Cfl(Bi-Fi)] 

= - [f j(B2-S)tf2(B3-Fj)+f3(B4-F2)] 


= + ^ f B 

2i 1 1 

= iV ^2 ®i 

= + ^ f B 

2i 3 1 


(B.9) 


The solution for the coupling and phase, mismatch are 
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62 J^r,2 «2 4 

K K ■ 12 L 1 " 16 ■ 10 " • • • 


(B. 10) 


(B. 11) 


where it is understood that \ = ± i y. The approximations to 1^ in- 

volves keeping only the first term in (B. 10) for r| « 1) and the first 
two terms in the series in (B. 11). Thus, for the sawtooth- wave. 




Ake 


K 


2 

_ A- 
12 


(0. 953) 


(B. 12) 


The approximation that uses only the first term in both (B. 10) and 
(B. 11) can be gotten directly from (3. A24, 3.A25) with ~ 5^ accuracy. 

At the third Bragg order all three of the periodicities con- 
sidered here contain the Fourier component f^. Hence, the major 
cross-coupling will be due to T] f ^ and the major self-coupling will be 
due to r| f j. The phase mismatch and coupling from (3. A24, 3.A25) 
are 

h. 

. K 


Ak £‘ 
K 


256 


(B. 13) 


X3 _ S 

K " 8 


(B. 14) 


The error is again on the order of 1^ except for the phase mismatch 
of the sawtooth- wave where the error ~ 5^. ^ 

The results of this appendix and comparisons to the sine wave 
are found in Tables B. 1 and B. 2. 
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Table B. 1 Normalized coupling coefficient Xn' ^ sine 
square- triangular- and sawtooth- wave periodicities. Fo 
the sawtooth- wave the coupling is x^ = ^ iX* 
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Table B. 2 Normalized phase mismatch 6j^/K for sine-^square 
triangular- and sawtooth-wave periodicities. 
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APPENDIX C 

REFLECTION AND TRANSMISSION 
COEFFICIENTS OF DFB FILTER 


Consider the following coupled equations 

Fj(z) - i 6p F,(z) = iXpqB,(z)"\ 

-Bj(z)-i6qBj(z) = iXqpFj(z) j 

which describe waves in a DFB filter. The primes denote differen- 
tiation with respect to the coordinate z. Differentiating (C. 1) and 
eliminating Bj(z) and Fj(z) respectively, produces 

HI I ■ 

V(b) = (C-2) 

where the subscripts and arguments have been dropped for simplicity 
and the primes denote differentiation with respect to the coordinate z. 


The solutions to (C. 2) are 

6-5 


F e 


where 


B e 


D 


6 -6 


r Dz , , -Dz 

= f j e + e 

u J- u -Dz 

= bj e + e 




(C.3) 
(C. 4) 

(C.5) 


Apply the boundary conditions 


F(i/2)= 1 


F(Z/Z)= T 


B(-i/2) = R 


B{Z/Z) = 0 


(C.6) 
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to find the following equations 


6 '-6 




-Di/2 f Di/2-| 2 

e + f^ e 


T - [fj e°^/2 


-[h 


• Di/2 


-Di/2-| i( ^2 )-^/2 

+ f 2 e Je 

Di/2"1 -i( 2~ 

+ 6^6 I e 


6 -6 


^ r Di/2^, -Di/21 

0=L^ie ■*‘^2® 


(C. 7) 


(C. 8) 


(C.9) 


(C. 10) 


Use (C. 7-8) to solve for f^ and use (C. 9-10) to solve for 

2 2 


[ ^(i 6i/2 ^ Di/2) _ ^ ^(-i 6 i/2 ± Di/2) "I 

Di TdI J 

^ e - e 

[- 6 in I DJ/2) 1 

h = ? |_ ^Dl _ ^-Dl J 


where 6 = (6 -6 )/2 

P q 


(C. 11) 


(C. 12) 


Two other equations for fj and bj can be found by substituting (C. 3,4) 

2 2 

into the coupled equations (C. 1). This produces the coupled equations 


(i 6 ± D)f, - i 6p £, = i Xpq bj 


(C. 13) 


-(i 6±D)bj - i-6^bj= i Xqpfj 
2 2 2 


(C. 14) 


Equations (C. 11-14) provide the necessary relations to solve for T 
and R, the transmission and reflection coefficients. After algebraic 
manipulations, the results are. 
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R 


T 


1 X 


32 . 


6 + 6 ^ 


D coth DJI - i(— 

i(6 -6 )Ji/2 
D e P ^ 

(6 +6 ) 

D cosh Dl - i — P-. ° sinh DZ 


(C. 15) 


(C. 16) 


No.te that the equations hold for all Bragg orders when the 
proper 6 and \ are used. Also the equations hold for coupling betv/een 
waves of different transverse modes, or in general for any coupled 
system where the phase mismatch and coupling might be different for 
the waves F{z) and B(z). 


V 
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APPENDIX D 

APPROXIIVIATIONS FOR DFB THRESHOLD AND SPECTRUM 


1. High-Gain Approximation 

Oscillation will occur in a DFB laser when the reflection or 
transmission coefficient becomes infinite. This implies the condition 


D^ ^ ' ^N 


(D. 1) 


where 


D 


N 


r 2 . 2 1 2 

= ^N ~ ^ 




^N ■ ^ ^N 


Equation (D. 1) can be written in exponential form 


D.,, - i 5 


'N 


N 




e = - 1 


(D.2) 


which is a complex equation. Under the high-gain approximation 
gj^ » Xj^ expanding the expression for Dj^ we find 


^N ^ ^ ^N ■ °N ^ ^ ^N ''' ®N 


(D.3) 


Either root of (D.3) when substituted into (D.2) yields 


(8n + i 


-(gN'*’ 


1 X 


= T 


N 


(D. 4) 


Taking the amplitude of (D. 4) we find the threshold gain condition 


for g 


N 


.,2.2 * ^®N^ 

^(8n ^N ^ " ^N ^N ® 


(D.5) 
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where the value is found from equating the phase of (D. 4), or, 
tan"^(Aj^/gj^) - Aj^i = (m + + phase (x) (m=0, ±1, ±Z, . . . ) (D. 6) 


This latter result produces multiple values for the spectrum or oscil- 
lation frequencie'^s which are called longitudinal modes. Since 
tan” ^(Aj^/gj^) « 1, the mode spectrum can be further approximated 
by 


Akyi” 

r 

K 


- (m+i)„ - (^^2 ) ^ 


2 ^t3 
T) N 


(D.7) 


where 


C„) •{ 


32(N‘'-1) 

0 index coupling (all N), gain coupling (N even) 
ir/Z gain coupling (N odd) 


for index or gain coupling and singly periodic media. 

The high-gain approximation yields the entire spectrum aind a 
transcendental equation for the threshold gain of each longitudinal 
mode. 


V 


2. Low-Gain Approximation 

It is more convenient to start with an alternative expression 
for index coupling or even- order gain coupling. Consider the solutions 


Fi =fi e^^ + f^ 


, D£ + b- e 
Bj = bj e 2 


-DA 


to the coupled equations 


Fj-i6NFi = ixNBi 
-Bj-i =iXiqFj 


(D.8) 

(D.9) 

(D.IO) 


where 
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Since there are no outside sources for DFB lasers, we may use the 


boxindary conditions 


B(- A/2) = 0 = F{X/2) 


This implies the relation 






When (D. 12) is substituted into the coupled equations (D. 10) along with 
expressions (D. 8, 9), the following equation can be derived. 


sinh = ± i D. 


Either of the exact equations (D. 1) or (D. 13) may be used for threshold 
calculations. 

Consider the low- gain approximation where gj^ « index 

coupling and even-order gain coupling. Expand the hyperbolic func- 


tion in (D. 13) to obtain 


Xn + 




where 




From the real part of the above equation. 


A 2 2 - ,,.2 


For large, oscillation takes place just outside the bandgap edges 
given by 


(D. 16) 


From the imaginary part of (D. 14) the gain condition becomes 




(D. 17) 


The low-gain approximation for odd-order gain coupling is 
most easily derived from equation (D. 1).’ First consider the neces- 
sary condition on the perturbation for no average gain (i. e. g = 0). 
An expansion of the hyperbolic function in (D. 1) produces 



(D. 18) 


Equating the real and imaginary parts produces the results 


= 0 (D. 19) 

2 

|Xn\ I =3 (D.20) 

Thus oscillation takes place at the bandgap center. The necessary 
length for a given perturbation, is the critical length 




(D.21) 


For gain perturbations or lengths less than those given by (D.20), no 
oscillation takes place. 

To find the threshold condition on the average gain, equation 
(D. 13) can be expanded to yield 








Keep all terms in to get 


3 ! 


± i 


(D.22) 


-172- 


and' 






0 


(D.23) 


(D. 24) 


Hence oscillation takes place at the bandgap center with the approxi- 
mate gain condition 


when i is large. 

Note that the low-gain approximation gives explicit values for 
the spectrum and thresholds of 'the average gain and perturbation for 
the lowest order (m=0) longitudinal mode only, 

23 

All previous results match those of Kogelnik and Shank at the 
first Bragg order. This appendix extends their results to all Bragg 


orders. 
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APPENDIX E 

ECW EQUATIONS FOR TRANSVERSELY PERIODIC MEDIA. 


Consider the wave equation 


+ 


dz 


dx 


+ k^e(x) 


E(x, z) = 0 


(E. 1) 


where e(x) = e (1 -+t]cosKz) for transverse singly periodic media. 
Assume a TE field of the form 

E(x,z) =y%(z)j_2„/N e'C<l-2"/N)PoX + ipx+k^o»] 

+ n(- 7 \ i[-(l-2n/N)PoX+Apx+kz z] [ {E. 2) 

+ (J,)S(z) 

for N even 
for N odd 

for N = 1 

for N even 
for N odd 

This accoxints for positive group velocity space harmonics involved in 
cross- and self-coupling of Fj(z) and Bj(z). A similar expression 
exists for negative group velocity space harmonics. 

The following set of N+3 coupled equations exist when (E. 2) 


where 


N"/2 


N-2 

2 


/ N-1 

"" 2 


Q- 


V. 

1 

0 


-1 


is substituted into (E. 1) for higher-order Bragg interactions. 


-174- 


tk^.-k,/ - (p„ + .p)"} F, + 2ik,^F/ = .i!|a,F,^^.Fj.3^ 

( J )tk^e-k4- 4p^]S = .( o>4^(F2/n+^;n) \ (E.3) 


{k . -k,^ - (-p^+.p)^ ] Bj +2ik,^ b; = . 


fl'^'-kzi-Kl«/N)Po+''P])En-2A^«ikz„B’j^2/N= »! 




where 


kz_ + ^k 


,P 

k 


+ Ap = NK/2 + Ap 

k + Ak 
o 


k^e = k^^ + P^ 

o o 

The primes .denote differentiation with respect to the coordinate z 
and the arguments with respect to z have been dropped for simplicity. 
The above equations are solved using the ECW assumptions and ap- 
proximations of chapter III. 

First, solve for F(z).j £/N 2/N* restilt is 

2 

similar to (3. A16-17) with the change T) -• T|/sin 0^ where 
1. 

sin9 = P /k e^. Hence, ... , 

° ° ° , ,>N+1' tl ,, 

-C T 1 N F (z) '■** V 2„ J 

' 1 ' . 2 sin 0Q 


F(z) 




‘-2/N ■ Pd-l/Njsin^e.^ JJ " (4„(„.n)/n2} ^ 


(E. 4) 
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B(z) 


1-2/N 


- Cj^ Ti N Bj(z)- 
~ 


+ 2«n_0o 


8(1-1/N)sin J|- (4 „(„.n/n2)2 


where 


fO 

■ { 1 


N = 1 
N ^ 2 


r (n-d/2 

f"'(n) for N odd 

n=l 


TT*fV)= ^ 


N/2 

■yy f(n) for N even 


n=l 

1 


for N = 1 


Second, trivially solve for ^(^)2+2/N ^^^^1+2/N* result is 

again identical to the previous resxilts of (3.A18-19) with the change 
2 

n -« n/sin 0 . 

' o 

_ C2,tiNF^(z) 

1+2/N 8(1+1/N)sin^e^ 

ri N B (z) 

®(^)i+2/N ~ 2 

L+d/SS 8{1 + 1/N)sin^e^ 

Substituting (E.4-7) into the equations for Fj(z) and Bj(z) in (E. 3) 
produces the ECW coupled equations for TP media. 


(E.6) 

(E.7) 


Fj(z) - i 6j^ Fj(z) = i Xjg Bj(z) 


Bj(z) + i Bj(z) = iXj^Fj(z) 


(E.8) 


where 


-176- 






2k 


(E.9) 


, ,.N+1 , 2^^N 

(-1) k e T| 


{sin^e Tf*[4n(n-N)/N^}^ 


(E.IO) 


The expressions can be further simplified by noting the equality 


k e^/kz = l/cos0 


2 2 2 “ 
k e = 0 + kz and the relations k e^/kz 

o "^o o o u o 

1 i 

and k e-^/0 = l/sin0 . The deviation from k and S represent the 

o "^o o . . o *^o 

changes in frequency and angle from the exact Bragg condition. 
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